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Abstract: Theories which have fundamental information destruction or decoher- 
ence are motivated by the black hole information paradox where one appears to 
have pure states evolving into mixed states. However such theories have either vi- 
olated conservation laws, or are highly non-local. Here, we show that the tension 
between conservation laws and locality can be circumvented by constructing a rela- 
tional theory of information destruction. In terms of conservation laws, we derive 
a generalisation of Noether's theorem for general theories, and show that symme- 
tries imply a strong restriction on the type of evolution permissable. With respect 
to locality, we distinguish violations of causality from the creation or destruction 
of space-like seperated correlations. We show that violations of causality need not 
occur in a relational framework, although one can have situations where correlations 
decay faster than one might otherwise expect or can be created over spatial distances. 
This creation or destruction of correlations cannot be used to signal superluminally, 
and thus no violation of causality occurs. We prove that theories with information 
destruction can be made time-symmetric, thus impossing no arrow of time. 
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1. Introduction 



All of our current theories of nature are unitary - an initially pure quantum state 
evolves into another pure state. Evolution is completely predictable and reversible. 
But it doesn't have to be this way - one can construct theories in which pure states 
evolve irreversibly into mixed states - information is destroyed, and predictability 
breaks down. If we follow the dictum: "That which is not forbidden is required", 
then we ought to seriously consider such non-unitary theories unless they are ruled 
out by other considerations. Indeed, that such evolution should be considered on this 
ground alone has long been advocated [[I], [2|], but it was the black hole information 
paradox || that led to serious consideration of non-unitary theories. Another moti- 
vating factor may be found in some interpretations of quantum measurement, since 
some interpretations involve non-unitarity in the form of fundamental or intrinsic 
forms of decoherence, or spontaneous wave-function collapse. Even in some versions 
of the many-worlds interpretation, the effective description of one of the branches is 
an effective non-unitary theory. 

The motivation coming from black hole physics is most easily seen when we create 
a black hole from an initially pure quantum state and then allow the black hole to 
evaporate away into thermal radiation Without radical changes to physics as we 
know it, the evaporation process must be non- unitary. Either Hawking's calculation 
holds, and the pure state evolves into a mixed state of thermal radiation, in which case 
the evolution must be non-unitary; or, Hawking's semi-classical calculation breaks 
down, and information does escape from the black hole through subtle correlations, 
in which case the evolution is again non-unitary. This is because if evolution takes 
a state outside of the light cone (from A to B for example), and our theory is 
relativistically invariant, then there exists a set of hypersurface [|5| in which the state 
has evolved from an initial copy at A, to two copies of the state, one at A and one 
at B. Such an evolution cannot be unitary — it violates the no-cloning theorem || 
which states that no unitary (or even linear) map can take any unknown state ip and 
produce two copies of it: tp ® t ^ ip ® ip. In the case of the black hole, one finds a 
space-like hypersurface which is well-away from the singularity, yet intersects almost 
all the outgoing Hawking radiation as well as the infalling matter which formed the 
black hole inside the apparent horizon. This hypersurface contains two copies of the 
state. Thus if information eventually escapes the black hole, the no-cloning theorem 
(and hence, unitarity and linearity), would be violated. We thus have the amusing 
situation that if no information escapes from the black hole, unitarity is violated, yet 
if information escapes from a black hole, unitarity once again appears to be violated. 

This is not necessarily a paradox - we are merely forced to conclude that funda- 
mentally, our laws of evolution might be non-unitary. The other possibilities, most 
notably, that the casual structure of space-time might get modified even in the semi- 
classical regime, is very radical, although there is some evidence for this in the form 
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of the AdS-CFT conjecture [[/]]. It might also be that all the information escapes 
at the very end, when the black hole has shrunken to the Planck scale, becoming a 
remnant with arbitrary high entropy and all the problems this entails H, M (but c.f. 



In this article we do not take a strong position on whether the laws of evolution 
are non-unitary - we simply stress that it is a possibility, and therefore the conse- 
quences ought to be explored much more fully. It is not enough to reject non-unitary 
evolution based on aesthetic grounds. All of the current candidates for a theory of 
quantum gravity, such as string-theory or loop quantum gravity, are unitary theories, 
while a serious discussion of non-unitary theories is strikingly absent. 

That a full exploration of fundamentally non-unitary evolution has not happened 
yet, can be traced in large part to a result by Banks, Peskin and Susskind (BPS) [fLl]] . 
There, the authors argued that non-unitary evolution must either be non-local, or it 
must violate conservation of momentum and energy. 

This tension between conservation laws and locality has not been overcome, and 
has posed a stumbling block to further exploration of non-unitary theories. While 



not denying that this tension existed, Unruh and Wald|12| argued that if the non- 
unitary part of the evolution was at sufficiently high energy, then the violations in 
momentum and energy conservation would also be at high energy, and hence, not 
observable in the lab. An attempt to implement such a theory in a natural way has 



been proposed simultaneous to the current work fllB] . Unruh and Wald also noted 
that if the evolution laws had a memory, one could restore the conservation laws. 

We will not adopt these approaches here, although we feel they have merit and 
deserve further study. Rather, we will argue that locality and conservation laws 
are not fundamentally incompatible. We will begin in Section |2| by discussing the 
objection of Banks, Peskin and Susskind, and the general form that non-unitary evo- 
lution must have. We note that their objection also applies to many interpretations 
of measurement in quantum theory. We then show that locality and conservation 
laws can exist in harmony provided one considers relational theories. For example, 
non-unitary theory which decoheres a particle into its position with respect to an ab- 
solute notion of space will violate momentum conservation, but one which decoheres 
two particle into relative positions will not violate momentum conservation, and can 
still be local. 

For pedagogical reasons, we first discuss momentum conservation and relational 
degrees of freedom in the context of quantum mechanics in Section ||] before moving 
to the full field theory case in Section |j. However, from a conceptual point of view, 
many readers will find Sections 0-|3| satisfactory for learning the central idea behind 
our construction. Next, we clarify different notions of locality. In particular, we 
note that one ought to distinguish theories which violate causality, and the far less 
problematic form of non-locality where correlations might be created or destroyed. 
In particular, we find that the relational theories we are interested in can lead to the 
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creation of short-lived correlations between two spatially separated regions due to 
indistinguishability of particles in relational theories. The existence of these correla- 
tions need not violate causality or allow super-luminal signalling. Additionally, one 
can have situations where the decay of correlations is slightly faster than what one 
might expect from our intuitions about non-unitary theories. Again, there does not 
appear to be any troubling consequences due to such effects. Demonstrating that 
one can have non-unitary theories which satisfy conservation laws and are causal is 



done in Section |5.1| Then in Section |5^2| we show how they may cause correlations 
to behave differently than in Hamiltonian evolution. We thus find that non-unitary 
theories are not non-local in the traditional sense, but rather have a benign and 
subtle form of non-locality which does not lead to violations of causality and appears 
to be a small effect. The theory respects a form of minimal Lorentz covariance, in 
that the equation of motion is invariant under the proper orthochronous Lorentz 
group. Lorentz invariant models in the context of non-unitary theories where first 



introduced in [ H . In Section |5| we prove that non- unitary theories must violate a 
stronger form of Lorentz invariance, in that the unitary violating term cannot be a 
Lorentz four-vector. 

At this point, it is worth mentioning that from a mathematical point of view, 
one can alway take a non-unitary theory, and make it unitary by considering an 
enlarged Hilbert space|l| where one adds an environment. However, in Section [5] it is 
apparent why the theories described here should be considered as fundamental. The 
environment that needs to be used is one which doesn't exchange any momentum, 
energy or other conserved quantities with the system. It is also in a highly non-local 
and artificial state, involving an infinite number of fields, and the interaction must 
also be engineered in a very particular way. So although one could consider it to 
be a physical environment which is invisible to ordinary interactions, this would be 
an unnatural view, since the environment interacts with the visible world in a way 
which is completely unlike any interaction which commonly exist, and incredibly 
complicated and engineered. It is thus far simpler in such a case, to consider the 
non-unitary theory to be fundamental. 

Having shown that causal non-unitary theories can conserve momentum, we next 
show how to construct theories which conserve any other conserved quantity. In con- 
trast to the other conservation laws, conservation of energy needs to be treated 
differently. This is for a number of reasons - in particular the Hamiltonian is no 
longer the generator of time translations, and as a result time-translation symmetry 
no longer implies as strong a restriction on the dynamics as do other symmetries. 
Nonetheless, we show that time-translation symmetry if applied to the entire uni- 
verse, does require conservation of the total energy of the universe. In Section [F| we 
discuss these restrictions, and construct a local theory which satisfies them. We also 
discuss the differences between Energy conservation/time-translation invariance and 
other conservation laws, and symmetries. The latter are treated in Section R 
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In Subsection 8.1 , we discuss the question of how rich are the set of observables 
which can describe relational quantities locally. This is relevant because these are 
the set of observables into which our theories can decohere a system into. We claim 
that this set, which is also the set of observables which commute with all conserved 
quantities, is a large set, and can thus lead to highly non-trivial decoherence. Another 
way to say this, is that the Hamiltonian, if it describes the full system and not just 
a part of it, is highly degenerate and there are thus many observables which can 
commute with it. 

One should however ask whether we ought to attempt to satisfy conservation 
laws. This is because Noether's theorem only applies to unitary theories. For non- 
unitary theories, it is not clear whether a symmetry implies a conservation law. We 
address this question in Section |9|, and derive a generalization of Noether's theorem 
for non- unitary theories. Surprisingly, one finds that symmetries do in fact, lead 
to strong constraints on the dynamics, although exact conservation laws may be 
relaxed in some circumstances. These results apply not only to theories which are 
fundamentally non-unitary, but also to effectively non-unitary theories which occur 
when a system interacts with an environment but still possesses some symmetry. We 
derive the resulting continuity equation and modified conservation laws for the stress 
energy tensor. When coupled to gravity, stronger constraints on the stress-energy 
tensor may also occur, and hence one might have stronger constraints on conservation 
laws. 

We then turn to the original motivation for considering such theories - and 
discuss how to couple such theories to gravity and in particular, black holes. We 
mention a few routes one might take in such an endeavor in Section TH . We also note 
that because non-unitary theories can create correlations without violating causality, 
they open up possible alternatives to inflation, for solving the horizon problem (the 
fact that the universe appears homogeneous, even though distant regions in space 
are causally disconnected). This is discussed briefly in Section pi]. 

Next, we turn to another motivation, commonly cited for non-unitary theories, 
namely, theories of stochastic collapse. Such theories attempt to explain the so-called 
"collapse of the wavefunction" by appealing to a stochastic mechanism. While we 
do not necessarily advocate such attempts, we do note that the relational theories 
considered here can solve some of the difficulties which have plagued stochastic col- 
lapse models. Namely, such models violated conservation laws and one could not 
construct a field theoretic version of such models. We show how to do so in Section 
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We conclude in Section [12] with some of the many open questions which remain. 
Throughout the paper, we try to keep the discussion as general as possible, but 
in Appendix [X] we discuss a number of particular model theories which may be of 
interest. We also show how correlations behave in some of these models in Appendix 
|R In Appendix we discuss a question which arose in the context of our discussion 
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on Lorentz invariance, namely, whether such non-unitary models have an arrow of 
time. It is generally claimed that this is the case. However, we show that one can 
construct a time-symmetric version of any non-unitary theory. Given the state of a 
system at some time t , predicting the future evolution of the state in a non-unitary 
theory is just as difficult as retrodicting its past, which is natural in any theory with 
randomness. 



2. The tension between locality and conservation laws 

The state of a system in both quantum mechanics and quantum field theory, is 
described by a density matrix p. Any new evolution laws must evolve the density 
matrix into another density matrix, and must therefore be trace preserving and 
completely positive. They must also be linear in the density matrix if the meaning 
of the density matrix is to be maintained, since a probabilistic mixture of two states 
ought to evolve into a probabilistic mixture of each of the evolved states. The 
evolution laws are thus completely-positive, trace-preserving (CPT) maps A. We 
will generally restrict ourselves to theories which do not possess a hidden memory 
- the evolution law at time t depends only on the state of the system at time t 
and not at other times (the evolution is said to be Markovian, or is generated by a 
semi-group) 

A(t x )A(t 2 ) =A(t l + t 2 ) (2-1) 

This feature can be motivated by the fact that if our theory is fundamental, there 
are no additional degrees of freedom which could be used to store the information 
about the system's past. However, one can have theories which are non-Markovian, 
due to a non-locally in time if one looks at very short-time scales. Discussion of such 
theories, as well as other relaxation of the conditions we impose, is beyond the scope 
of the present work. 



It can be shown that the most general Markovian map for bounded operators] 16 



17fl is given by 

^ = -i[H, p] - l - J2 Hj{M}M iP + pM^Mi - 2M lP M}) (2.2) 

»j 

with 7y a positive matrix and H is the usual Hamiltonian. Since 7^ is positive, 
one can always consider it to be diagonal. This is referred to as the Kossakowski- 
Lindblad master equation, or usually just Lindblad equation fl8|]. The operators 
are usually referred to as Lindblad operators, and we will adopt this convention also. 
We will sometimes refer to the Mj as observables when they are Hermitian, because, 
as we will see, the resulting evolution given by Equation ( |2.2| ) is very similar to the 
evolution of a system after measurement of the Mj. The term on the right hand side 
in addition to the Hamiltonian commutator is sometimes called the dissipater T>(p). 
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The sum over operators could be replaced by a continuous spectrum of operators, 
in which case the sum would be replaced by an integral. In the bulk of this paper 
we will consider the discrete case, since extension to the continuous case is usually 
st r aight f or war d . 

Two examples are of particular interest. The first is when the Mi are projectors 
Pi = and 7^ = 7^ in which case the Lindblad equation decoheres the state p 

in the \i) basis, with off-diagonal matrix elements decaying exponentially fast 

to zero at a rate exp — ft. This is pure decoherence, and is much like a measurement 
in the \i) basis (although without collapse into a particular |i}). A second example, 
is to take only one operator M, and have it be a self-adjoint operator (i.e. an 
observable). This will also cause decoherence in the eigenbasis of the operator M, 
with a decay exponent for matrix elements \j)(k\ given by (mj —rrik) 2 , with rrij being 
the eigenvalue of M in eigenstate In the case of Hermitian Lindblad operators 
Mi after diagonalization of 7^, we can rewrite the Lindblad equation as 

i 

Mathematically, one can think of such maps as being the result of a unitary 
evolution acting on the system and a hidden environment. However, we shall see 
that for the theories considered here, it would not make any physical sense to think 
of it in this way. 

Banks, Peskin and Susskind (BPS) rediscovered the master equation (|2.2|) , and 
argued that its evolution, if local, will lead to a violation of both energy and momen- 
tum conservation. We will not reproduce their result here, but rather give intuitive 
arguments which contains the core ideas. A particular example of this Lindblad 
evolution would be the proposal of Hawking [H] which BPS recast as 

§ = -iW,p] - 2^4 / dx[F^F, v {x\ [F™F Ta (x),p]} (2.4) 

where we have simply taken the Lindblad operators F^ v (x) to be the electromagnetic 
tensor. Such Lindblad operators are local, and the evolution law leads to massive 
production of energy and momentum, even if the initial state is the vacuum. It 
is not difficult to see why. The operators F^ix) do not commute with the total 
momentum P of the field. The evolution corresponds to the system interacting with 
a local environment, and each interaction gives the particles a kick of momentum. 

We can see this clearly in the Heisenberg representation - the evolution of oper- 
ators G which are conserved by the usual Hamiltonian evolution, will instead evolve 
under the master equation as 

dC 1 

~dJ = -^HjiMjMjG + GM\M 3 - 2M-GMj) (2.5) 
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If we demand that G is exactly conserved, then the Mi ought to commute with it. 
However, when G is the momentum P of the field, only highly non-local operators 
Mi can commute with P. Thus, BPS argued that non-local operators Mj will result 
in non-local dynamics. 

There is a connection here with quantum measurements. Closely related to the 
evolution of Equation fl2.4p - not in the context of field theory but within quantum 
mechanics - is to consider Lindblad operators which are projectors onto the position 
x, P x = \x)(x\ leading to evolution 

^ = -i[H, p] - 7 / dxP xP P x (2.6) 

Again, one can think of the particle as interacting with a local environment at each 
point x, and the momentum of the particle will be disturbed. Such a master equation 
leads to decoherence which one can think of as very similar to a measurement. The 
evolution given by Eq. ( ^.6[ ) decoheres a particle into position eigenstates. Although 
there is no collapse, it is as if the position of the particle is being measured. A 
measurement of position disturbs the momentum of the particle, since 

[X,P}^0 . (2.7) 

In much the same way as Hawking's evolution Equation (|2.4| ) interacts locally and 
thus disturbs the momentum, so a position measuring device interacts locally with 
particles and generally disturbs their momentum (c.f. fT"9|j). 

There is thus something disturbing about the conclusions of BPS. Although they 
are considering evolution laws, the exact same considerations would seem to apply 
to measurements. Depending on one's interpretation of measurement, one may be 
forced to conclude that momentum is not conserved. However, we locally measure 
the position of particles all the time, and we don't believe we are violating momentum 
conservation. But if P is changed when we measure X, how is it that we can measure 
the position of a particle without violating momentum conservation? 

The answer is that we don't actually measure X of the system with respect to 
some absolute reference frame, but rather we measure X$ of the system relative to 
the position of some measuring apparatus Xr which acts as a reference system f20|1 . 
We can thus measure the relative position X$ — Xr, without disturbing the total 
momentum P$ + Pr of the system plus reference, since 

[X s -X R ,P s + P R ] = (2.8) 

This suggests that we should consider Lindblad operators which are relational. We 
will now see that relational Lindblad operators can allow one to have local interac- 
tions which conserve momentum. Before moving to the modified evolution laws in 
the case of quantum field theory, it will be instructive to first look at the problem in 
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the context of ordinary quantum mechanics. After using this as a toy model we will 
proceed to evolution in the context of quantum field theory, which will allow us to 
examine locality issues in greater detail. 



3. Relational evolution within quantum mechanics: local dy- 
namics with momentum conservation 

As an example, we wish to produce a theory which decoheres particles into position 
eigenstates. Normally, this would be given by Lindblad operators of the form P x = 
\x)(x\ resulting in Equation ( |2.6| ) above. Since this violates momentum conservation 
let us add a second particle (a reference particle) and consider instead the Lindblad 
operator 

Q= / dx\x 1 ){x 1 \ <g> \x 2 )(x 2 \ (3.1) 



and evolution of the form of Equation fl2.2|), but with non-relational Lindblad opera- 



tors Mi replaced by the relational operator Q. In analogy with operations sometimes 
performed in the context of quantum information theory [21], we will say that Q 



is created from by group averaging. In quantum information theory, this is 

sometimes called twirling. We will henceforth use Q, to denote a Lindblad operator 
which is relational, while general Lindblad operators will be denoted by Mj. We will 
also have need for local Lindblad operators in field theory. These we will denote by 
Li(x). 

Using Q in place of M leads to evolution of the form 

d JL = -i[ H , p ] - i 7 (Q p + P Q - 2QpQ) (3.2) 

Q is the relevant operator for coincidence measurements, i.e. it corresponds to the 
observable that two particles are at the same position. One can think of the first 
particle as being the detector, and the second particle as being the system (or visa 
versa). Q then measures whether the detector "clicks". By using this relational 
operator Q as a Lindblad operator in Equation |2.2j we get a theory which only 
decoheres two particles when they are in the same place. At all other times, the 
operator Q does nothing, and one has ordinary Hamiltonian evolution. This evolution 
is thus local. On the other hand, total momentum P tot = P\ + P 2 is conserved. This 
is easy to see by rewriting Q = J dxe~ lPtotX \x ) (x \ ® \x ) (x \e lPtotX which clearly 
commutes with the translation operator e~ lPtoty (Q is actually invariant under it's 
action due to the bounds of the integral being from minus infinity to positive infinity), 
and hence [Q, P tot ) = 0. 



Let us now modify the Hawking model of Equation fl2.4|) to ensure momentum 



conservation. To do so in the context of a quantum mechanical reference system, we 
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add a non-relativistic object. Rather than taking F^ v (x) as the Lindblad operators, 
we take the relational operator 

Q = J dxF^F tMV {x)®\x){x\ (3.3) 

and evolution of the form 

f = -i[H, P ] - ^Kr<U [Q Tr7 Qr„ P ]} (3.4) 

Because of the integration over all space, it is again easy to see that the Q^ u commute 
with the total momentum of the field and reference system. It is also local - the 
difference between this evolution and the initial one is that decoherence only occurs 
at the location where the reference system is, and is zero everywhere else. If the 
reference system is at a particular location such that the total state of the system is 
initially 

P = P field ® \xq){xq\ (3.5) 

and the particle is heavy enough that the wavefuntion doesn't spread fast, then 
Equation |3]4| gives 

df jf L = -i[H,p\ - £^[F^F, u (x ), [F™F T(7 (x ),p fteld }} (3.6) 



which looks exactly like Equation ( [2. 4|) for the field degrees of freedom. 

Note that the above evolution law may not be too far from what we are aiming 
for. It produces decoherence at the place where our reference object is. Although a 
black-hole is not a particle, one might imagine that it acts very much like a quantum 
particle for the purposes of these considerations and that something like the above 
evolution law might be what we are seeking with being approximately the 

position of the singularity. 

It should also be noted, that although total momentum is conserved, the momen- 
tum of each system need not be. For the relational Lindblad operator of Equation 
( |3.1|) , when both particles are coincident, and hence undergoing decoherence, they 
will get a momentum kick in opposite directions. This can be lessened either by 
making 7 sufficiently small, or by smearing the projection operators, for example, 
instead of M x = \x)(x\ <8> |^)(#|, we can take it to be a smeared projection 

M x = J dyg(x-y)\y)(y\®\x)(x\ (3.7) 

with g(x — y) either Gaussian exp [—A (re — y) 2 ]or some other peaked function around 
x. This means that there will be violations of causality on the scale of the width of 
g(x — y) (this is the same degree of causality violations which appear in the original 
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spontaneous collapse theories which we discuss in Section 11). The momentum ex- 



change between individual particles can be lessened more easily in a field theory, a 
point we will return to in Section |j. However, on a fundamental level, one should not 
be concerned about momentum exchange (it of course happens in unitary theories 
as well), except to the extent which it may violate experimental constraints. 

Within quantum mechanics, we can perform group averaging on any local oper- 
ators, either smeared or not, the idea is just to build in some redundancy with extra 
particles so one has non-trivial Lindblad operators. We now move to a discussion of 
similar dynamical models as the one given above but in the context of quantum field 
theory. Once we have shown that one can conserve momentum while still preserv- 
ing locality in field theories, we will move to the more complicated case of energy 
conservation. 



4. Relational evolution for fields: momentum conservation 

Consider a raw local Lindblad operators Li(x) equivalent to the local operators im- 
plied by Hawking's requirements. 

Because these are local in absolute space x, they cause large violations of energy 
and momentum conservation. To conserve momentum, we proceed as in the single 
particle case - we introduce a second field, and consider relational Lindblad operators 

Qi,(t) = J dxLf\x,t)®Lf\x,t) (4.1) 

and evolution of the form 

^ = -i[H, P ] - ~ 'M/hji' + pQljQij - iQaPQlj) (4.2) 

ij 

The operator Qij(t) clearly commutes with the total momentum just as the Q of 
Equation |3.1| did. We will see in the next section that the evolution is local. We 
will take 7^ to be diagonal (we could have just as easily considered the non-diagonal 
case) . 

It is also not hard to see that in the non-relativistic limit, evolution of the form 
of Equation (|4.2|) can yield the evolution of Equation (|3.2|) described in the context 
of ordinary quantum mechanics. We show this in Appendix |A[ 

We can consider more general forms of dynamics than that which would be 
given by Equation ( |4.1| ). For example, one could use only one field and make it 
"self-referential" , as is done in one of the models in Appendix [A]. Or we can consider 
any local operator L^x) which can act on as many fields as we wish, or be a product 
operator or not, and then perform group- averaging it to make a relational Lindblad 
operator Qi(t) 

Qi(t) = [ dxLi{x,t) (4.3) 
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and evolution of the form of Equation ( f4.2|) Again, this clearly conserves momentum, 
for the same reason as the two previous examples (i.e. Qi(t) is invariant under the 
action of e lPx ). 

Note that although the total momentum is conserved, there will be momentum 
exchange between the two fields under such models. While this is not an issue 
from a fundamental perspective, one would like to have some control over it. We 
can do this by making the Lindblad operators less projective. One procedure is as 
follows: for a local Lindblad operator L(x), let us first smear it over a small test 
function in space which has support on a very small point in space. This is the usual 
smearing over a tiny region which just keeps the theory finite, and we will generally 
assume in this paper that all field operators are suitable smeared. We then take the 
spectral decomposition of the smeared L(x), and we call the projectors Pi(x) onto the 
eigenvalues I of L(x). We can now consider another set of local Lindblad operators 

Ui(x) = J dl 9l (l)Pi(x) (4.4) 

with gi(l) a highly peaked function around some value of I. This new set can be as 
degenerate as we wish (for example, if ^(/) is a constant function, then IT (a;) will 
just be proportional to the identity. It will thus not disturb the individual systems 
at all. By choosing g^l) appropriately, we can control the momentum exchange by 
making the operator IT(x) more or less degenerate via tuning g^(l) to be more or 
less peaked. We can then use IT (2) or some function of it, in our group-averaged 
Lindblad operators. E.g. for two fields 

Qij = j dxTlf\x)Tlf\x) (4.5) 

We now investigate the theory's locality properties. We will find that the theory 
respects causality if the Qi satisfy some simple conditions. 



5. Locality of the field equations 

Although BPS did not specify what they meant by locality, it will prove useful to 
consider two different notions. The most important one is that the equations of 
motion ought to be casual, i.e. they should not allow super-luminal signalling. In 
this sense, we will see that the evolution of Equations ( fL2l - f4~3l) is local as long as 



it satisfies certain conditions - for example, that the Lindblad operators Li(x,t) 
are local, and that they are either Hermitian or come in pairs with their Hermitian 
conjugate. In Section |6] we will see that not only does the theory satisfies causality, 
but also a form of Lorentz invariance. 

However, as noted in the Introduction, the causal evolution of the fields is not the 
only aspect of locality. There is also the question of how correlations evolve. As in 
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! 1 l| we need to consider how spatially separated operators behave under non-unitary 
evolution. We will see that correlations can be created between space-like separated 
regions, however since these correlations cannot be used to signal between two dif- 
ferent regions at speeds faster than light, there creation doesn't violate causality and 
doesn't overturn any cherished physical principles. 

The effects of the creation of distant correlations are also experimentally benign 
- the only consequence is that in special circumstances, correlations decay slightly 
faster than one might expect, or can be created. Some constraints from causality on 
how correlations can change in the context of EPR experiments has been explored 



in 
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We will first in Subsection |5J] show that causality can be preserved by construct- 
ing a model which is derived from a fully Lorentz invariant field theory by tracing 
out some degrees of freedom. Since the original theory is causal (and even Lorentz 
invariant), and the Lindblad equation is derived merely by ignoring some of the fields, 
the Lindblad equation clearly respects causality; ignoring degrees of freedom cannot 
result in super-luminal signalling. The construction works for any Lindblad equa- 
tion whose operators Qi are integrals of local field operators. The field theory and 
associated environment are very artificial (and apparently, necessarily so), and thus 
should not be thought of as an actual fundamental theory from which the Lindblad 
equation is merely an effective description. 

Next, we explore, In Subsection |5.2| how correlations evolve. The issue of Lorentz 



invariance is taken up in Section ^ 

5.1 Causal and relational evolution 

Causality requires that for local field operators 

[A(x,t),B(x',t')\ = (5.1) 

when (x—y) 2 < 0, and the evolution of A(x, t) and B(x', t') are given by the Lindblad 
equation. We shall now show that relational evolution of the form of ([4.2| - f4.3|) does 
not necessarily imply a violation of relativistic causality. More precisely, we show 
that this evolution is causal to at least second order in the Born approximation. To 
this end, we consider an explicit construction of a Lindblad model starting from a full 
model of a system and a (fictitious) environment which together are fully relativistic 
and therefore causal. 

Consider then a number of fields <fii(x) which are coupled to some environment 
of relativistic fields ip a {x) via some Lorentz invariant interaction. The full reduced 
dynamics, obtained by tracing out the environment must therefor be causal. 

Our construction takes an interaction between the fields given by 

Hj = / L(x)A{4j Q {x))dx, (5.2) 
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where for simplicity we only take one local operator L(x) which acts on the (pi, and 
with A(ip a ) denoting some local operator built of the environment fields {ipa}- One 
can consider many different Li by just including more interaction terms. 

The equation of motion for the full density matrix in the interaction picture is 
given by 

dp(t) 



dt 

with an integral form of 



+[H t (t),p(t)] (5.3) 



p(t) = p(0)-i [ dsiHjis^pis)] 
Jo 



(5.4) 



which can be inserted into Equation ( p.3| ). Tracing out the environment degrees of 
freedom in which the ip a (x) live, and associated with the density functional pE yields 
the equation of motion for p$, the density functional of the fields <f>i(x) 



dp s (t) 



f dstrEiHrit),^)^^)]] (5.5) 
Jo 



dt 

where we have assumed that 

tr E [Hj(t),p(0)}=0 (5.6) 

We will see that we can make the interaction weak, and in such a way that 

tr E H lP (t) « tr E Hj (p s (t) ® p E ) (5.7) 

Note that this condition is often written as 

p(t) « p s (t) ® p E (5.8) 

but strictly speaking there will be correlations between the system and environment; 
it is just that the interaction term is not sensitive to them. An example of this is 
an atom coupled to a field where the correlations between the atom and field (in the 
form of emitted photons) are carried off to infinity and don't interact again with the 
atom. 

Inserting the form of the interaction Hamiltonian, Equation ( |5.2| ) and the approx- 
imation Equation ( |5.7[ ) into the equation of motion for the reduced density matrix, 
Equation ( |5.5| ) then becomes 

dp s {t) 



dt 



ds \ dxdx'(A(x,t)A(x',s))[L(x,t),[L(x',s),p s (s)]]. (5.9) 



This is in fact the second order Born approximation and accounts only for the lowest 
order corrections to unitarity. Interestingly however, in the particular case that the 
state of the system and environment is Gaussian the above form can be shown to 
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be exact [23[|. We then take t to infinity to remove explicit dependence on the initial 



poo p 

I dsJdxdx f (A(x,t)A(x,t-s))[L(x,t),[L(x',t-s),p s (s)]]. (5.10) 



time. This approximation is valid as long as Equation (|5.7| ) holds approximately (i.e. 
the system is Markovian enough). This also will be shown to be justified by suitable 
choice of the interaction and environment. Rewriting the above integral then gives 

dps{t) __ 
dt 

We must now depart from the standard microscopic derivation of the Lindblad 
equation. It is easy to see that consistency with momentum conservation imposes 
requirements on the environment field correlations. If we demand translation invari- 
ance (i.e. Lindblad operators of the form Equation (|4.3| ), then the environment needs 
to have constant spatial correlations. We also need fast decay of temporal correla- 
tions to reduce Equation ( |5.10| ) to the Lindblad form of Equation ( |2.3|) . Namely, we 
require 

(A(x,t)A{x',t- s)) =75(t-s) . (5.11) 
Equation (|5.10| ) then yields the traditional form of a Lindblad equation with Hermi- 



tian Lindblad operators 

with the momentum conserving Lindblad operator being identified as 

Q= L{x)dx . (5.13) 



Clearly any choice of a local system operator L(x) will also conserve momentum. 
The question is then can a Lorentz invariant theory manifest correlations such 
as in Equation ( |5.11|) ? Under familiar situations, such as the vacuum or thermal 



states, the field correlations decay both in space and time, thus the environment 
here cannot be taken as a familiar field state. However, we can design an interaction 
and environment with such properties by considering an environment with an infinite 
number of fields ip a at hand. We can then tailor an interaction in such a way that at 
each instant in time, the field interacts effectively with only one environment field: 

A(x,t) ~ J da8(t- a)A a (x,t). (5.14) 

with A a (x, t) a function of ip a only. As a result, for t ^ s, the temporal correlations 
disappear and we have (A(x, t)A(x', s)) = (A(x,t))(A(x',s)). We can always choose 
the field and interaction terms such that (A t (x,t)) = while (A t (x,t) 2 ) ^ so that 
Equation (|5.11| ) is satisfied. 



We further need to ensure that there are constant spatial correlations in the 
environment. We can choose the environment state to be in a superposition of states 
which are constant over space 

\if> a ) = / dtpg a (tfj a )Il x \tfj a ) x , (5.15) 
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where \ip a )x denotes the field state at the location x. This is sufficient to enable 
Equation (|5.11|) to be satisfied. For example, one can take A a (x,t) to be ij) a (x, t) 
and g a (ip a ) chosen symmetric around if) a = 0, so that (A(x,t)) = while (A t (x,t) 2 ) 
is a constant and independent of x. 

The free environment Hamiltonian is irrelevant for us and can be set to zero for 
simplicity. Since the system interacts with different degrees of freedom at each time, 
the temporal correlations vanish and the Markovian approximations we have made 
can be justified. Note however that the Lindblad equation only arises up to second 
order in perturbation theory. It is possible that we might not get an exact Lindblad 
equation from a causal model. However, our model very closely approximates Lind- 
blad evolution, thus the relational Lindblad evolution is causal, or at least closely 
approximated by a theory which is. Up to this degree of approximation, we see that 
our relational theory of Equations ( p~2| - [4.3| ) can be derived from a fully causal and 
Lorentz invariant theory, and is thus causal. 

Note that the construction above works for Hermitian operators Qi, but one 
can construct an interaction Hamiltonian which produces non-Hermitian Qi, as long 
as the Qi come in pairs with their Hermitian conjugate. This is achieved via the 
interaction 

Hj= [ (L{x)A\i)i(x)) + L\x)A{i)i(x)))dx (5.16) 



with A(ipi(x)) a bosonic operator. We then require, instead of Equation ( J5.ll ), that 
(A(x, t)A(x', t-s)) = (A f (x, t)A f (x', t-s)) = Q 

(A(x, t)A\x', t-s)) = (A\x, t)A{x', t-s)) = j6(t - s) . (5.17) 

This can be done for example, if A(x, t) is a complex field operator, and then charge 
conservation will require that terms like (A(x, t)A(x', t — s)) vanish. We can then use 
a similar procedure to get the Lindblad equation 

% = -\i(Q ] Qp + pQ f Q ~ ZQpQ ] + QQ ] p + pQQ ] - 2Q^ P Q) (5.18) 

where again 

Q = dxL{x) . (5.19) 



Equation ( |5.18j ) is of the form of a Lindblad equation, with two Lindblad operators 
Q and (JL Thus, using a relativistic environment to derive causality, we get either 
Hermitian Lindblad operators, as in Equation (|5.12|) , or non-Hermitian Lindblad 
operators that come in pairs as in Equation (|5.18|) . We will see in Section |9] that 
these are pairs of annihilation and creation operators. In Equation (|5.18|) , we have 



only one pair of Lindblad operators, but by having many terms in the interaction 
Hamiltonian, one can have a Lindblad equation of the form 

^ = -\n{Q\QiP + pQ ] Q l - WiPQl + Q 1 Q\p + pQ%Q\ - 2Qt p g i ) (5.20) 
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If we now try to find a more general condition which will guarantee that a 
given Lindblad equation is causal is complicated by the difficulty in solving the 
evolution equation for all times. However, we can more easily explore the question 
of whether instantaneous signalling is allowed. This is just the requirement that the 
evolution laws of a local operator should not depend on distant observables. I.e. 
dA ^ = C(A(x)) needs to be a local operator. The conditions for this to be true are 
explored in Appendix [D]. These conditions clearly includes the ones used here, but 
may be more general. 

5.2 Evolution of correlations 

Let us now look at the second notion of non-locality, namely, how correlations may 
evolve. Here we show that the effect of the non-unitary evolution is fairly benign. 
Consider two spatially separated operators A(x) and B(y), potentially acting on 
many fields. To control divergences it should be understood that we would actually 
integrate these operators over some test function highly peaked around x and y, but 
we shall omit this smearing for ease of notation We want to check what happens to 
correlations between these two local observables. If the evolution does not change 
the nature of correlations between distant locations, then we would expect that a 
product observable would evolve product. I.e. 

j t [A(x)B(y)] = j t [A(x)}B(y) + A(x)±[B(y)) (5.21) 

Local Hamiltonian evolution H = J dx7i(x), clearly satisfies the above condition. 

Let us now look at how correlations evolve in the non-unitary theory given by 
the momentum conserving relational model of Equations (f4.2[ - [4.3|) , and taking again 



L(x) to be a local operator. We ignore the Hamiltonian term since that evolution is 
purely local. 

£[A{x)B{y)] = -\Y,liAQ)QiA{x)B{y)+A{x)B{y)Q)Q i -2Q)A{x)B{y)Q j ) 



f dzdzXL](z)L i (z')A(x)B(y)+A(x)B(y)L]{z)L i {z^ 



2 

■t 



-2L}(z)A(x)B(y)L 1 (z')) 
jA{x)B{y) + A{x)jB{y) + V(A(x)B(y)) (5.22) 



where V is the terms which violates the product evolution form of Equation ( 5.21 
and is given by: 



V{A{x),B(y)) 



Y^HjJ dz[Li(z),A(x)] J dz[Lj(z),B(y)] (5.23) 
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while the local terms are 
j t A(x) = -i^ 7w (L}(x)L i (x)>l(x)+A( i r)L}(x)L l (x)-2Lt(x)A(x)L i (x)) (5.24) 

and similarly for B(y). Because is a local operator, J dz[Li(z) , A(x)} and 

J dz[L_j(;s), B(y)} in Equation ( |5.23| ) will be a local operator at x and y respectively. 
Note at this point that A(x) evolves locally and does not depend on any observables 
or evolution at y. This is enough to guarantee no superluminal signaling at an 
instant, a point which is discussed in greater detail in Appendix [TJ. 
Putting it all together, we have 

j t [A(x)B(y)] = j t [A(x)}B(y)+A(x)j t [B(y)]-Y,^ f dz[L t (z), A(x)} f dz[Lj(z),B 

ij 

(5.25) 

It is noteworthy, that if we allow distinguishable particles, such as the models 
discussed in Section [|, there is no violation of the product rule (Equation (|5.21| ) ) 



for the evolution of observables. This can be seen as follows. Take as the Lindblad 
operators, those of the kind given by Equation (|3.3|) 



Qi = J \z)(z\Li(z)dz (5.26) 

and imagine we have one distinguishable particle. Then the local observables which 
evolve via the disipater are only those which contain operators which act on this 
particle such as \x)(x\A(x) and \y)(y\B(y). However, it makes no sense to speak of 
the evolution of \x) (x\A(x)\y) (y\B(y) - this term is zero unless x = y. 

This has potentially deep consequences in a fully relational theory. In a relational 
theory of distinguishable particles, one can describe observables at two points in space 
by relating them to two distinguishable particles, one at each of the two points. 
Furthermore, in a non-relational theory, one can, in effect, distinguish two identical 
particles by distinguishing where they are. Measurements on an electron on earth 
commute with measurements made on an electron on the moon. However, in a 
relational theory of indistinguishable particles, one is bound to have a violation of 
Equation ( |5.21| ) which merely reflects the fact that in a relational theory, we can't 
use absolute space to distinguish two points, and if the field values are identical at 
those two points, then there is nothing which can distinguish those two points. As a 
result, the evolution laws must act on both those two points without differentiating 
between the two of the. 

Another interesting consequence of Equation (|5.25 ) can be seen if one takes the 



limit that x — > y and A = B. We then see that under non-unitary evolution, powers 
of operators do not necessarily evolve as the power of the time evolved operator, i.e. 
the relation 

d/dtA n (x,t) = nA n ~\x,t)d/dtA(x,t) (5.27) 
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which holds for unitary evolution, need not hold here. We do not know of any 
dramatic consequences of the failure of this relation, but it might be worthwhile to 
explore such effects. 

We thus see that although the evolution laws do not violate causality, they do 
not satisfy Equation ([5.21 ) I.e. correlations between two regions might be created 
or destroyed. This will be shown more explicitly in Appendix [B] when we consider 
particular models. Such an effect cannot necessarily be used to signal superluminally. 
Thus, although the long-range creation of correlations might make us uncomfortable, 
the fact that it does not necessarily lead to violations of causality requires us to take 
such a possibility seriously. 

For the purpose of looking at the evolution of long-range correlations, the time 
evolution of the mutual information between two spatial operators is a good figure of 
merit, i.e. we are interested in I(A(x,t) : B(y,t)) where I(A : B) = S(A) + S(B) — 
S(AB) and S(A) is the von Neumann entropy of a quantum state pa, S(A) = 
—trpA log pa- Analyzing such a quantity for a particular model is beyond the scope 
of the present article, but would be of interest. 



6. Lorentz invariance and non-unitary evolution 

Thus far, we have not concerned ourselves with the question of Lorentz invariance, 
although we saw in Section || that causality was preserved. This was because we can 
view non-unitary evolution as coming from a fully Lorentz-invariant theory, with an 
environment which we trace out. The question we now address, is to what extent 
the new theory is Lorentz invariant - it may be that the traced out environment 
provides a preferred frame which breaks Lorentz invariance. 

There are various forms of Lorentz invariance one might demand: the first is 
that the right hand side of the Lindblad equation ought to transform like d/dt, so 
that both sides of the equation transform in the same way under a Lorentz boost. 
This has been referred to as a minimal Lorentz invariance requirement. It has been 
claimed |24| that even this minimal requirement is impossible to satisfy. However, 



while certainly difficult, one can adjust the construction of Lindblad operators and 
coupling constant in order to get an equation which is covariant in this first sense. 
This is done, for example in [T^, |14], [25j . 



In the context of relational theories, Lorentz invariance arises quite naturally, 



at least in terms of the minimal requirement demanded in p4 |. This will become 
apparent in Section [7|, when we see that Lindblad operators are Lorentz scalers if 
they conserve energy and momentum, and are of the form 

Q = dxoltL(x,t) (6.1) 
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This allows one to construct a minimally transforming evolution law 

^ = -i[H, p] - \ J] 7 ii(0]0iP + pQtQ, - 2Q\pQ j ) (6.2) 

where we take 7^ to transform as 4-. I.e. 7 become dynamical, as it is in any model 
where the non-unitarity is induced by tracing out an environment in a unitary theory. 
In a derivation based on tracing out the environment, 7 is related to environmental 
degrees of freedom. 

This ensures that both the left-hand side and right hand side transform under a 
Lorentz transformation in the same way. One can think of it as a redefinition of the 
time derivative with 

d d „ . , „ 



with T>q just the usual dissipater term taken from Equation ( |2.2|) and transforming 
in the same way as d/dr. Note that it doesn't matter whether we view this as a 
redefinition of d/dt, or d/dr. 

The requirement of this minimal Lorentz invariance appears to restrict the Lind- 
blad operators to those which are of the form of Equation (|6.1|) , and (as we will see in 
Section [?], thus imposes energy conservation). However, other forms are possible as 
well, at least if one relaxes the locality requirements. However, it certainly appears 
to give strong weight in favor of energy conservation. 



This minimal form of Lorentz invariance, introduced in |L4j, uses the Lindblad 



equation for the evolution with respect to a time-like vector a, i.e. d a , transforms 
like a Lindblad equation, while, for a space-like vector, b, one has the usual relation 

d b p=-i[P b ,p) (6.4) 

with P b the momentum in the direction b. For proper orthochronous Lorentz trans- 
formations, time-like vectors a remain time-like, while space-like vectors b remain 
space-like, and thus, this difference between spatial and temporal translations is pre- 
served. There is thus a distinction between time translation (which is not unitary) 
and spatial translations which are given by unitary operators, and are in fact con- 
strained to be unitary transformations. One might wonder whether the term T> Q (-) 



in Equation (|6.3|) can transform as a full four-vector. 

We now prove that it cannot, provided that the usual definitions of energy and 
momenta are unchanged. In the Heisenberg representation the equation of motion 
of a scalar field 4>(x) is then 

d Q( j>{x) = -i[<t>{x),H}+D Q { ( t>} (6.5) 

The invariance means that the equation must have the same form 

d Q $(x) = -i$(x),B\ + DS(Z)] (6-6) 
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in another reference frame with x = Ax. 

For contradiction, we now assume that Dq indeed transforms as a time compo- 
nent of a four vector D^, hence 

D-S(x)]=<Do{<P}+A 1 D 1 {<j ) } (6.7) 

where the boost is for concreteness taken along the "1" direction. 

In addition we also have 4>{x) = <f>(x) and H = Aq/J+AJPi and Bq = Aq<9o+Aq<9i. 
Substituting, we then get the necessary condition for invariance 

d 1 (^ = -i[<f>,P 1 ]+D 1 [ < f>] (6.8) 

We observe that the first two terms are equal due to our assumption that Pi is a 
momentum component of a four vector P^. I.e. d\<f> = —i[<f>,Pi] is a kinematic 
identity, which follows from the definition of P % as 

P* = / dxn^l (6.9) 
J dxi 

and the canonical commutation relations (in contrast to [(f), H] which has no such 
identity associated with it). Hence we obtain the requirement 

Di[0]=O (6.10) 

I.e. Do cannot be a component of a four vector, and we only have the minimal form 
of Lorentz invariance described above. This ends the proof. 

This is not particularly troubling, since as demonstrated in Section [| it doesn't 
lead to violations of causality. It is as if there is a preferred reference frame which 
distinguishes time from space. 

We comment that it is possible to choose a nontrivial Dq in such a way that 
Do[4>] — vanishes identically (e.g. Q = J <j)dx) but is non-zero in the general case 
Z^ol 71 ""] 7^ 0. In this case <p satisfies the ordinary relativistic equation. Nevertheless 
non-invariant corrections will then be present in the equations of motion for higher 
powers of (ft or its conjugate field it. 

For completeness, we note that applying transformations of Equation ( |6.3| ) to 
the Klein-Gordon equation gives for the time-like component 

(□ + ^ 2 )<P{x) = -V\V {cj>)) - 2V°(d <j ) ) (6.11) 

and the Dirac equation would become 

i^d^ -m^ = 7°£>oW (6.12) 

with the 7^ being Dirac's gamma matrices rather than the coupling constants of the 
Lindblad equation. 
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7. Energy conservation 

As we have noted, energy conservation is in a slightly different category to momentum 
conservation. We will show that it can be conserved exactly. However, it may also 
be that it only needs to be approximately conserved. There are several reasons for 
this - the first is that a priori, energy conservation is in a slightly different category 
to momentum conservation: since the generator of time-translations is now given 
through the Lindblad equation rather than by the Hamiltonian, Noether's theorem 
no longer applies - time-translation symmetry does not imply energy conservation. 
We will see in Section ^ that Noether's theorem is also modified in the case of 
momentum conservation, but the change is far less severe. The lack of a Noether's 
theorem is both a blessing and a curse - on the one hand, it implies that we only 
need to conserve energy approximately and only to the extent that experiments place 
constraints on energy conservation. On the other hand, without Noether's theorem, 
it is not immediately clear what prevents us from seeing massive changes in energy. 
As was noted in Section [| Lorentz invariance may provide some answer to this 
question (we will see why, further in this section). Additionally, constraints may 
arise due to locality, or due to coupling the theory to gravity, both of which demand 
energy conservation. We will discuss the former in Section |9| and the latter in Section 
10. 

For the moment, let us consider constructing an energy conserving Lindblad 
equation. Unlike other conserved quantities, this will involve two steps. The first, 
will be to construct group-averaged Lindblad operators as we might do for any other 
conservation law. In order to ensure locality, we will need a second step, which is to 
impose time-translation invariance. We proceed with the first step: we construct a 
relational Lindblad operator from any other Lindblad operator via, 

Q ( = /*e-<V. (7.1) 

and use the notation Qi to henceforth denote time averaging. 

Now, also as before, it is useful to have additional degrees of freedom which 
act as a reference frame. In this case, we need a system which behaves as a clock 
(although the distinction as to what is a clock needn't be artificial). To find the clock 
observable, we find a decomposition of H into Il r and H such that H = U T + H , 
with II T conjugate (or perhaps approximately conjugate) to some observable r, and 
such that r commutes with H Q . In that case, r will correspond to the physical time, 
as can be verified by it's equation of motion 

f = -i[H,r] 
« 1 (7.2) 

t may correspond to some global observable such as the expansion parameter of the 
universe; it may correspond to many local degrees of freedom which taken together 
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constitute a clock; it may correspond to a single field. For the moment we just let it 
be arbitrary, and potentially some collective degree of freedom. We then take 

Qi = j dte- iHt Qi <g> \0}(0\e iHt . (7.3) 

where |0) is the initial state |r = 0) of the clock and Qi is our original Lindblad 
operator (e.g. the momentum conserving relational operators of the previous section). 
It is not hard to see that Qi = Qi{r) i.e. it is equal to the original Lindblad operator 
at the time r as measured by the physical clock. 



Qi = j dte- iHot Qie iHat (g)\t}(t\ 



dtQi(t) ® \t)(t\ 
= Qi(r) (7.4) 

where the last line follows since Qi(r)\t) = Qi(t)\t). 

Any Lindblad operator given by Equation ( |7.1|) commutes with the total Hamil- 
tonian H, since it is invariant under the action 

e- iHe Qe iHe = Q (7.5) 

and therefore conserves energy when used in the Lindblad equation. However, such 
a quantity is rarely local in space. 

This is partly because of the explicit t dependence (in the form of a derivative) 
on the left-hand side of the original Lindblad equation. Additionally, although a 
field operator ip(x,t) commutes with ip(x',t) for x ^ x', the symmetric statement 
is not true: [ip(x,t),<p(x,t')] ^ for t ^ t'. This is crucial, since our derivation in 
Subsection | of the locality of the field operators relied crucially on the fact that the 
Hamiltonian was local, or, to put it another way, that 

dzLi(z),<l>(x)] = Dt(x) (7.6) 

with D^(x) some local function. We would not obtain a similar expression if we 
were integrating over time (c.f. pEfl ). The Qi are not local in terms of the absolute 
background space. 

One can think of this problem in another way - Because the left hand side of 
the Lindblad equation is a derivative with respect to absolute time t, it is as if the 
Lindblad operator is being measured (or decohering) at some parameter time t which 
is typically not equal to r . As a result, the operator Qi{r) is generally smeared over 
all space. As an example, consider the operator 5(X(t))S(r(t)) - it is completely 
smeared over all space except at the time t = 0. If 5(X(t))5(r(t)) appears as a 
Lindblad operator in Equation (|2.2|) then it is as if the operator is being measured at 
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an arbitrary time t. What we really want is that this operator should be measured 
at t = r not at some other time t. We will now see that it is possible to have the 
former situation. 

Namely, we should now impose time-translation invariance. The entire density 
matrix of the universe cannot depend on the external time t i.e. this parameter time 
is physically meaningless. As a result, we must have 



This is similar to the Wheeler-deWitt equation p7| for gravitational systems except 
that the total energy of the universe is not constrained to be zero, but rather, some 
constant. Since the density matrix is stationary with respect to this external time, 
the total energy of the universe is constant. As a result, Equation [T7| is enough to 
guarantee that the absolute vacuum state, for example, is stable. 
It implies 

-i[H, p}-\j^ HjiQjQiP + p^A ~ 2 Q*pQ]) = ( 7 - 8 ) 

Decomposing H in terms of the clock Hamiltonian, and the rest, gives 

z[n r , p] = -i[H a , p\-\Y1 ~u(Q'/Ip + pQjQi - iQiPQ]) (7.9) 

If we consider the density matrix as a function of the physical observable T, i.e. 
p(T,T'), then we can define the total derivative in terms of the c- number r 



dp(T,T>) _ dp(T,T) , dp(T,T') 



dr dT dT 

Since z[II T ,p] = dp/dr, we have 



(7.10) 



d £ = -i[H , p]~IJ2 y^QjQiP + pQiQt ~ ^pQ\) ( 7 - n ) 



i.i 



which is precisely the usual Lindblad equation but with the external, non-physical 
parameter time t replaced with the physical observable time r . In terms of the 
external t, the energy is conserved 

dH 

- = (7.12) 

and as long as we take the Lindblad operators to be the time-averaged Q of Equation 
(|7.3|) , we will also conserve energy with respect to the physical time, ie. dH/dr = 0. 

In terms of locality, the Lindblad Equation ( [7.11Q is just as local from the per- 
spective of physical clock time r as the momentum conserving Lindblad equation 
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H) was in terms of t. In Equation (4.2), we have dp/dt on the left hand side and 



Lindblad operators Q(t) on the right hand side, while in Equation ( [7.1 1|) we have 
dp/dr on the left hand side and the same Lindblad operators, except as functions 
Q{t) of r. The locality properties are thus the same. 

This proposal that we make is not without it's difficulties, however, these dif- 
ficulties are also present in any unitary theory as well. Namely, if we believe that 
quantum mechanics is fully consistent and applies to the universe as a whole, then 



Equation (|7. 7j) will hold, and one will need to find an internal observable to act as a 
clock. Likewise, if we believe that Hawking radiation can be explained by a unitary 
theory of quantum gravity, then the Wheeler deWitt Equation akin to (|7.7|) holds. 
In other words if we are trying to decide whether our fundamental evolution laws 
are unitary or non-unitary, then both options suffer from the same problems when 
it comes to energy conservation and time. 

In a unitary theory one encounters various issues of time and energy conservation, 
usually studied in the context of attempts to apply canonical quantization to gravity. 
In particular, although it will be possible to find an observable r which is at least 
approximately conjugate to some II T , if it is not perfectly conjugate, then there is 



some (arguably small) probability that the clock will run backwards [28]. There is 
however, a lot of freedom in how we decompose H into H Q and I1 T . Not only can 
we choose it so that r acts as a good clock, but we can also choose it so that the 
energy exchange between the clock and system is as small as possible. This provides 
an interesting new criteria for decomposing closed systems, for example, in quantum 
cosmology. 

Additionally, it appears that the state p does not evolve in time t. This is 
sometimes referred to as the "frozen time" formalism. In fact, this is not a good 
description of what is happening. Just as we group-averaged Q in Equation (|7.3| ), 
we could have instead group-averaged the density matrix to get 

p= [ dte- lHt p®\0)(0\e iHt (7.13) 



so that given any p and physical clock observable, we get a state p which satisfies 

dt 

There is nothing frozen about p - it simply represents a state where we are ignorant 
of the absolute value of t i.e. the initial parameter time t. However, a hypothetical 
external observer, who knows what the absolute time t is, would see the evolution 
as e~ l ^*~*°)p (g) |0)(0|e^*~*°). We should thus think of the state of the universe as 
actually being in one of the many states in the integrand of Equation ( 7.13j ), and 



evolving in time t as p(t) = e~ iH ^~ to ^p ® |0) {Q\e lH ^~ to \ It is not frozen at all. We 
just don't happen to know the starting time to and therefore, which of the p(t) in 
the integrand happens to describe our universe. 
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In conclusion, we have seen that we can construct a Lindblad equation which 
conserves energy and is still local. In order to do so, we not only needed Lindblad 
operators which commuted with the Hamiltonian, but also we needed the absolute 
time t on the left hand side of the Lindblad equation to be replaced by a physical 
time r. The theory has to be completely relational. 



8. General considerations for ensuring conservation laws 

With the exception of time-translation, the procedure for creating Lindblad oper- 
ators which respect any particular symmetry is straightforward. One considers lo- 
cally gauge invariant operators Li(x) (which we call the raw operators). For each 
symmetry that one wishes to respect, one can add an additional field which acts 
as a reference frame (although no artificial distinction need be made between the 
system and reference frame). We will see that one can also have self-referential 
models where a field acts as a reference for itself. One then considers joint opera- 
tors on all the fields La' 2 ' (x). Often, these will be chosen to be of product form 
Lj*j£' 3 "'\x) = Li(x)Lj\x)Lip(x)..., but this is not necessary. One then considers 
the group generated by all the generators G of the symmetries, and we average over 
this group 

Qa = J dgU{g)L^U\g) (8.1) 

with dg taken to be the uniform measure (i.e. Haar) over the group, and U(g) the 
unitary representation of the group. Usually U(g) is obtained by exponentiation 
of infinitesimal generators U(g) = exp—iGg, but this is not always the case, for 
example, with diffeomorphisms. 

One can also consider what we call coherent group-averaging. I.e. operators of 
the form 

Q a = J dgdg'U(g)L^ 3 -W(g') (8.2) 

we discuss these in slightly more detail in Appendix |A|. 

Note that as in the cases considered here, we are often averaging over a non- 
compact group, and thus the operator does not formally converge. For the purposes 
considered here, such operators can still be considered in the Lindblad equation. 
A review of this and other technical issues can be found in the review of in 
the context of the quantization of constrained systems. A review of general group 
averaging in the context of work on reference frames and quantum information can 
be found in p0|. 



Note that one might want to consider variations of this such as considering 

(1 2 3 ) (1) (2) (3) 

operators L\^ " (x) = P> (x)Pj (x + h)Pk \ x + ^)--- as long as the l v are small 
enough that the non-locality of the above operator will be unobservable. There are 
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several models of this form which we have studied and we discuss a few of these in 
Section [A]. 

In general, there appears to be a large set of observables that we can use as Lind- 
blad operators, and many local observables which upon integration give relational 
observables. It seems that one can include all powers of a local observable along with 
conjugate observable, e.g. 4> n (x)7r m (x), and then perform group averaging on them. 
It might be useful to study such models in more detail. 

8.1 Richness of relational observables 

The question of what physical observables can be described by integrals of local field 
operators is important in order to be convinced that one can get non-trivial dynamics 
from the models we have proposed. This will become even more important when we 
discuss stochastic collapse models, since there, one is attempting to describe all the 
observables we measure in the laboratory. 

If one considers Lindblad operators of the form of Equation ( |S.1| ), it is not imme- 
diately clear how rich their structure is i.e. are there enough of them, that they lead 
to non-trivial decoherence. If the only Lindblad operators where the total energy, 
total angular momentum, total electric charge etc., then we would only be able to 
decohere into observables which were conserved quantities. This would not be very 
interesting. The question is whether there are many other observables besides the 
conserved charges which are of the form of Equation ( |3.1| ) 

Since these Lindblad operators can also be understood as relational observables, 
one can ask whether relational observables fully describe the world around us. Cer- 
tainly one expects that relational observables do describe what we observe, since 
these are the only physically relevant ones - they are the ones we measure in the 
lab. Indeed, when it comes to conservation laws which are associated with a Gauss's 
law (such as charge and energy), we can only measure operators which commute 
with the conserved quantities |nj. Thus, the relational ones (which commute with 
these conserved quantities), being the ones we measure in a lab, must be sufficiently 
rich. However, here, we are also demanding that the observables be local. Again, 
local relational observables would appear to be sufficiently rich, since in the lab, 
we are always measuring local quantities. Given that in the lab we presumably are 
measuring local observables which must commute with conserved quantities, it is 
reasonable to believe that Lindblad operators which correspond to those observables 
are sufficiently rich. 

However, it is not completely obvious that such a set could fully describe the 
types of observations we make. Take for example a particle located at x and a second 
one at y. The relational description would be limited to saying that there are two 
particles which are a distance \x — y\ apart. This relational observable, which can be 
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written as 



Q— I dz\x + z)(x + z\®\y + z)(y + z 



(8.3) 



is not local, and thus cannot be used in our Lindblad equation. Lindblad operators 
of the form 



might be permissible, but only if I is sufficiently small. The theory is non-local on 
that distance scale. 

Now, if one has a pre-existing reference frame, i.e. a field <f>(x) over all space 
which has different values at all positions x, then certainly it is simple to use such 
a field to fully describe all additional fields in relation to it. For example, as a toy 
theory one might be able to use a suitably smeared local projector P<f,(x) onto values 
of the smeared field 4>(x) and the relational observables 



If 4> is such that it uniquely distinguishes x, then this acts like a good reference frame 
and the observables Qi^ will be able to describe any observables with respect to this 
reference frame. This includes observables of the form of Equation ( |S.3|) . However, 
without a pre-existing reference frame, it is not clear whether all observables can be 
described in local and relational terms. Perhaps it is too much to demand that a 
relational model be descriptive without a pre-existing reference frame. 

Indeed, in a lab setting, we build up our knowledge of the reference frame over 
a period of time. A spark in a spark chamber (a local coincidence measurement) is 
equivalent to a position measurement only because the positions of the wires have 
been predetermined by making many measurements over a duration of time (e.g. 
sending photons to bounce off various parts of the device, and timing their return). 
It is thus certainly feasible that local relational observables, over a period of time, 
will decohere a system in a non-trivial manner, and eventually, completely decohere 
a system into states which, from a physical point of view, are the states we observe. 

In order to pursue such questions further, the Hamiltonian dynamics is likely 
to be as important as the Lindblad operators. For example, if there is an effective 
potential a(x)a(y)V(x — y) between two spins, then the probability of the two spins 
being parallel or anti-parallel will be a function of the distance between them. Lind- 
blad operators which decohere the spins are in a sense decohering different distances, 
particularly if there are many of them. Such models, and the general question of 
richness of local relational observables is one which deserves future study. 

Finally, we often try to find an observable which commute with all symmetries, 
and whose commutator with a local observable is a local observable - namely, when 
trying to find the Hamiltonian of a system. The Hamiltonian is an example of an 
operator which commutes with conserved quantities, however, renormalizability and 




(8.4) 




(8.5) 
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boundedness place more constraints on a Hamiltonian than it does on an observable. 
For a Hamiltonian, one is forced to throw away odd-powered terms like 3 which are 
unbounded from below, but such terms can be included in an observable. Likewise, 
for terms which would lead to non-renormalizable theories, it appears that one can 
include these in an observable. However, one needs to ensure that including such 
terms as Lindblad operators would not result in a non-renormalizable theory, or 
lead to other problems. Understanding this is an interesting question which we feel 
deserves further study. 

It is worth noting that for a physical system, if one includes perturbations due to 
other systems, then the Hamiltonian is usually non-degenerate. The energy levels, 
unless protected from external perturbations by a symmetry, will be split by the 
presence of the external fields. Thus, decohering a system in the basis of the Hamil- 
tonian is as decohering as one can get, and the Hamiltonian encodes the full richness 
of the theory. Knowing the energy is a complete observation of the system. However 
for closed systems, the Hamiltonian is generically expected to be highly degenerate. 
It is this fact which allows us to make all sorts of measurements, which still commute 
with the Hamiltonian. 



9. Symmetries and conservation laws 

Noether's theorem states that for a system which evolves unitarily, there is an equiv- 
alence between symmetries of the evolution laws, and the conservation of some quan- 
tity. When the evolution is no longer unitary, the connection between symmetry and 
conservation laws may break down. Most symmetries are probably regarded as being 
more fundamental than the associated conservation laws - if the laws of physics were 
discovered to be different on Earth than on the moon, we would be unlikely to con- 
clude that the laws of physics changed just because we moved to some new location 
- rather we would more likely posit that the presence of the moon or some unseen 
matter was causing objects to behave differently. We therefor ought to demand of 
our non-unitary theory that it be invariant under our cherished symmetries (such as 
spatial and temporal translations), but we may be less bothered by the fact that it 
may violate a conservation law. 

However, this presents a potential problem for non-unitary theories |32| since 
quantities like momentum appears to be conserved experimentally, and it would be 
difficult to explain such a fact without Noether's theorem. Without a connection 
between symmetries and conservation laws, there appears no reason to demand even 
approximate conservation laws. The near exact conservation laws we experience 
would therefore appear to be an unexplainable accident in a theory with non-unitary 
evolution. We will return to this point towards the end of Section |9.1| , 

However, although the connection between symmetry and conservation laws may 
break down for a general map which takes mixed states to mixed states, it need not be 
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the case. One has to look at the structure of the map to make such a determination. 
We will therefore re-examine this connection for the general evolution law given by 
Equation (EO) 



C(p) = i± = -i[H,p]-±Y, 7ij(MjMip + pMjMi - 2M lP M}) . 



'j 



We will find that symmetries do still place a constraint on the type of evolution 
which is permissible, although the connection with conservation laws is modified. We 
also make the distinction between symmetries which are associated with a generator 
of a unitary transformation (such as spatial translations), and a symmetry which 
is not. This is important, because in the theories considered here, time-translation 
is no longer generated by the Hamiltonian, but is instead given by a non-unitary 
operator. 

Let us consider symmetry under an infinitesimal transformation implemented by 
e~ lGe with G a Hermitian operator which generates the transformation. Symmetry 
of the evolution laws implies 

C(p) = e iGe £(e~ iGe pe iGe )e- iGe (9.1) 

I.e. the equations of motion "commute" with translation under G - if we translate 
under G and then evolve the system, and then translate back, the evolution should 
not be different. Stated another way, if we implement a transformation on the evolu- 
tion laws, and also transform to state, then if there is a symmetry, the transformed 
state will satisfy the transformed evolution laws. 

When applied to the master equation of ( |2T2l in the (potentially degenerate) 
eigenbasis \g) of G, invariance under the symmetry yields 

-iH-^j.MjMA \g)( g '\ + \g)(g'\ (iH-~'£ li M}MA +Y^ li M i \g)(g'\M} 
e9 e lGe i^-iH - ij^JWjMiJ \g)(g'\ + e^'\g)(g'\ \iH - i£ 7i JWjM<J e" lGe 



t„-iG?e 



+ e^'-^J2ne tGe M l \g}(g'\Mje 

i 

This can only be satisfied if 

[±iH-~J2liM}M i ,G}=0 (9.3) 

i 

and additionally 

e iG Mi = Mie iG - iA > (9.4) 



(9.2) 
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which ensures cancellation of the final term on both sides of Equation ( |9.2|) when 
used in conjunction with the Hermitian conjugate of Equation 



M\e~ iG = e~ iG+i ^Ml . (9.5) 



Here the Aj are real constants. Equations fl9.3H9.4Tj can be written as the two condi- 
tions 

[G,M i ] = -A i M i (9.6) 

and 

[G,H]=0 (9.7) 



with condition ( |9.6|) implying 

[G,M}] = A i M} (9.8) 
Another way to write the condition of Equation (|9.6|) is as 

YtiM^WWi = Y, m ^9X{g)n\9 + ^i)(g' + A,| (9.9) 

i i 

Since the effect of the Lindblad equation is completely characterized by its effect on 
a complete basis, we can always take each Mi to individually be proportional to a 
raising or lowering operator of G 

M}\g) = mt(g)\g + A i ) (9.10) 

and hence we get Equation (|9.6|) . This also implies that M/Mi commutes with the 



generator of the symmetry. 

[G,M/M,] = (9.11) 

thus M\Mi has the same eigenbasis as the number operator. Note again that \g) 
might be highly degenerate, and the values m*(g) can depend on other degrees of 
freedom - we have just dropped them for convenience. 

The Lindblad equation in the Heisenberg representation Equation ( |2.5|) , together 



with conditions ( p.6H9.7| ) gives 



f = J>Mt[G,M ( ] 

i 

Y. -,A,,U/.U, (9.12) 

i 

i.e. the rate of change of G is given by an average over the number operator. We 
may then write the solution to Equation ( |9.12| ) as 

G(t) = G(0) - (^ 7 ,AM/M 4 )t (9.13) 
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Here, G(t) could grow without bound, although it grows at a linear rate, rather 
than exponentially fast. From our discussion on causality, we found that the re- 
lational Lindblad equation was causal when the Mj were Hermitian, or when they 
came in pairs with their Hermitian conjugate. If we take these two cases, then we 
have that either 

G(t) = (9.14) 

as was the case when the fictitious environment was made up of scalar fields. Or, 
in the case of complex environment fields, then we have that for each Mj which is 
a raising operator by the amount Aj, there is a lowering operator Mj by the same 
amount. This gave a Lindblad evolution such as that of Equation ( [5.18| ). This implies 
that in Equation ( |9.12p , we must also put in the conjugate terms, which contribute 
a — AiMiMj to the sum (since they are lowering operators), giving 

G(t)=Y i7i A i [M},M i ] (9.15) 

i 

Comparing Equation ( |9.14| ) with Equation ( |9.15| ), one might want to rule out 
the latter possibility on physical grounds. This might be reasonable, given that non- 
Hermitian Lindblad operators require a different type of environment and interaction 
Hamiltonian (fictitious or otherwise). And that such an evolution would eventually 
lead to infinite production of the quantity G (although the time it would take may 
be very long). 

The only exception to our generalized Noether's theorem is for energy which is 
not generated by the Hamiltonian. Rather Equation ( |2.2| ) is itself responsible for 
time-translation, and is clearly time-translation invariant provided H and the Mj 
have no explicit time-dependence. Thus the only conservation law which requires 
additional explanation is that of energy conservation. We saw however, in Section |7| 
that imposing time translation invariance on the entire state of the universe led to 
energy conservation. We will also see in the next section that coupling to gravity 
appears to impose additional constraints. And finally in Section |6] we saw that a 
weaker form of Lorentz invariance also appears to imply energy conservation, at 
least for the type of models we considered. 

9.1 The modified continuity equation 

We now wish to see how the continuity equation is modified when evolution is non- 
unitary. If we have a generator of a symmetry which expands as G = J g(x)dx then 
under ordinary unitary evolution, we would have 

dg(x) 



dt 



i[H,g(x)\ . (9.16) 



Noether's theorem then gives that dG/dt = 0, which implies that J must be 
equal to a surface term which vanishes as the surface is taken to be at infinity. By 
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Gauss's law, this implies 

i[H,g(x)] = V i f i (9.17) 

i.e. the integrand of J dx[H, g(x)} can be written as a total divergence. For unitary 
evolution, this implies the usual continuity equation 

d„r(x) = (9.18) 

where one defines f°(x) = g(x). 

For non-unitary evolution, Equation ( |9.16| ) gets modified, and we have 



i[H,g(x)) -\Y,^ M l M ^)+9{x)M}M J -2M}g{x)M 3 ) . (9.19) 



dgjx) 
dt " l "" r/J 2 



From the symmetry conditions of Equation ( |9.7| ) derived in the previous section, we 
still have [H, G] = and thus Equation ( |9.17| ) still holds and 

i[H,g(x)] = V i f i (9.20) 

For local evolution, we can define 

^0) = ~J^l ij (MjM j g{x) + g(x)M}M j -2M}g(x)M j ). (9.21) 

From the results of the previous section, we know there are only two cases - 
either the M* commute with G in which case, once again, k(x) can be written as a 
total divergence, and one still has a continuity equation of the standard form. This 
will be the case if Mj are Hermitian. On the other hand, without this requirement, 
our generalization of Noether's theorem implies that the M l might also be raising or 
lowering operators of G. In such a case the continuity equation gets modified to 

d,r(x) = k(x) (9.22) 

and applying Equation ( |9.12| ) to Equation ( p.21[ ) we have that 

dxk(x) = ^^AigLlMi (9.23) 



As before, two cases are of special interest due to their locality properties. The 
first is when the Mj are Hermitian, in which case we recover the ordinary continuity 
equation. The second case is when Mj comes in pairs with it's Hermitian conjugate, 
in which case Equation (|9.15|) gives that 

J dxk(x) =Y^ 7i A i [MlM i ] (9.24) 
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as is implied by Equation Q9.15Q . However, in this case, in addition to the Mj coming 
in pairs with M i , we also required that the M« be integrals of local fields in order for 
the equations of motion to be causal. 

M i = J dxLi(x) (9.25) 

Because L(x) is local, the integral of the commutator [Ll, L(y)] is some local function, 
call it Ni(x)5(x — y). We then have that 

[M}, = [ dxdy[L\(x),U{y)} 

dxNi(x) (9.26) 

Putting this into Equation ( 9.24 ), we thus find that k(x) is a local quantity 

k{x) =J2^Mx) (9.27) 

i 

For the generators of spatial-temporal translations, we have a modified conser- 
vation equation 

T^. v = fc" (9.28) 

with 

k( x ) =J2lAiNi(x) (9.29) 

i 

for the /i which are the spatial components. I.e. since our version of Noether's theo- 
rem gives no restriction on energy conservation (since H is no longer the generator of 
translations), we can impose no restrictions on k l . However, the spatial component, 
k^, which we shall denote by k(x) must satisfy Equation (|9.23| ) with the possibility 



that all the Af are zero (the case of Hermitian Mi, in which case one has ordinary 
conservation of T^ v '. u . Note that here, the Aj is now a vector, since it is determined 
by Equation ( |9.6|) which depends on the symmetry generator G. As an example, Ni 
might be the number operator for momentum vector p. In this case, the Aj which 
may imply that coupling to gravity imposes the condition k(x) = 0. We discuss 
this briefly in Section [10]. It may be that curvature can not couple to matter in the 
usual way. There would also be a modification to the Gauss's law since energy need 
not be conserved in an asymptotically flat universe. Since the gravitational Gauss's 
law comes from Einstein's equation, it would necessarily get modified due to the 
modification of Equation (|9.28| ). With regard to the Bianchi identity, a potential 



resolution may be found in the model discussed in Section [5[ There, one enlarges 
the Hilbert space to consider a unitary theory on the original system and an envi- 
ronment. In such a model, the violation of stress energy in the non-unitary theory 
could be understood as an exchange of energy and momentum with the (perhaps 
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fictitious) environment. One could then have full energy- momentum conservation on 
the extended Hilbert space, and apply Einstein's equations to the stress-energy due 
to both the system's matter, and also the environment's. We leave such problems 
for further study. However, it does raise the issue of how we couple our theory to 
gravity, a problem which we very briefly address in the following section. 

Note that in general one expects that the modification to the continuity equa- 
tions may be difficult to observe. This is because the are total raising or lowering 
operators and would not have a big effect locally. I.e. local functions will be super- 
positions over many eigenstates of G, and so acting a lowering or raising operator on 
such a state would not change the state locally, although the global state will change. 

The modification of the continuity equation will lead to a modification of the 
Ward-Takahashi identity. Rather than use the ordinary continuity equation <9 M / M (x) = 
in the identity, we would use d^f^ix) = k(x) which leads to 



for T a functional of the fields and S e the infinitesimal gauge transformation. We 
hope to explore the effect of this modification further. 

10. Coupling to gravity 

Despite the fact that one of our original motivations for this study was information 
destruction in black-holes, we have thus far only considered non-unitary evolution in 
the context of flat-space, without any coupling to gravity. In all likelihood, coupling 
such a theory to black-holes will require a quantum theory of gravity. Nonetheless, 
we can construct a few plausible toy models in order to understand how such a theory 
might couple. Here me suggest a few different possibilities. We will not explore them 
in any detail, we merely mention some possibilities for future research directions. 

(1) Coupling to the singularity: If we want the modification to ordinary 
physics to only occur at the Planck scale, or in the presence of black holes, than 
one is likely to have information destruction at the singularity of the black hole. 
A potential toy model for such a coupling could be something along the lines of 
Equations (|3.3| - |3.4j) with \x)(x\ projecting onto the "position" of the singularity - 
where by position, we understand it to be in terms of some perturbative expansion 
around flat space-time. 




(9.30) 




(10.1) 




(10.2) 
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While this is only a very rough toy model, it has some of the features one 
might expect from black holes which destroy information. This type of coupling is 
particularly tempting in light of the fact that the quantum black hole presumably 
has a huge number of microscopic degrees of freedom owing to it's large entropy 
This makes its plausible that it corresponds to some sort of distinguishable object. 
As a result, we no longer have even the mild non-locality discussed in Section |5], as 
discussed following Equation ([5.261). 



(2) Coupling to the curvature: 

There are many ways one can create Lindblad operators which couple to the 
curavture or the metric in some way, and we mention a few toy models. For example 

Q = J dx^f=gR{x)L{x) (10.3) 

where L(x) is any local field operator and R(x) is the space time curvature or some 
other function of the metric E.g. one could consider vector components of the form 

Qfiu = J dx^gR^ u L{x) (10.4) 



or 



Q = j dxV^gR^(x) (10.5) 

or any such combinations, and functions which need not be in product form. 

Such Lindblad operators would not conserve momentum or energy, but we do not 
expect conservation if we were to treat gravity as an external field. One would have 
to take into account the back-reaction on the gravitational field, thus such models 
should only be considered as an effective description of the matter degrees of freedom 
in a fixed background. 

Note that these models appear to have the required effect of producing deco- 
herence near black holes, while producing little decoherence in flat-space. This is 
because the curvature terms goes to infinity as we approach the singularity, and the 
decoherence terms will dominate, even for very small coupling constants 7. 

Other relational observables in the context of effective gravity has been discussed 
3"3"| and in the context of deSitter space times in [33]. Such observables can be 
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used in the relational Lindblad equation. For example, we can take 

Qr = J dx^f T (R(x)) (10.6) 

where f T (R(x)) is a highly peaked function around R = r and can potentially serve 
as a clock. 

(3) Black holes as microscopes to high energy physics: One can consider 
a very simple idea, where we don't directly couple gravity to the master equation, but 
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instead use the fact that a black-hole redshifts energy near the horizon. If the dissi- 
pater term T>(p) only acts at high energy, then in ordinary laboratory experiments, 
one would find evolution very close to unitary. On the other hand, modes which 
are radiating from a black hole were once at a very high energy since they originate 
from close to the horizon. They will therefore have been subject to the dissipater 
term T>(p) causing information loss. Whether such a scenario is compatible with the 
usual picture of black-hole evaporation is unclear, but the idea is perhaps attractive 
because of it's simplicity. 

(4) Information Destruction: We have seen that it is possible to construct 
relational theories with little constraints other than locality and conservation laws. 
This evolution may be constrained further, by going back to one of our original 
motivations for considering non-unitary evolution. The goal was to use it to destroy 
information in the context of black-hole evaporation. In order to remain consistent 
with the causal structure of the black hole space time, no information can escape the 
black hole until such time as the black-hole is of Planck size, at which point the causal 
structure may break down. This means that virtually all the information which goes 
into the black hole must be destroyed. We therefore demand of our evolution, that 
it be strong enough to completely destroy information. 

What we mean by this is that for p ou t,in the density functional of the fields inside 
and outside the black hole, we have 



| \Iout ® Knpout,in ~ Pout ® T^lW - 6 ( 10 - 7 ) 

log a 

for all p ut,in and some very small e, i.e. we may want to allow for some small 
deviations from this the map being completely randomizing. Note that this condi- 
tion of completely randomizing is stronger than demanding that the map be merely 
randomizing, i.e. that 

hnPin = 7^ (10.8) 

log a 

We should also exclude information stored in global charges. This is because 
quantities like the mass, electric charge and angular momentum are measurable out- 
side the black hole and are conserved during evaporation. Amusingly, it is easy to 
use relational Lindblad operators to destroy all information except the information 
that is stored in global charges, because the Lindblad operators commute with these 
global charges. It is only the relational information one wishes to destroy. 

In the case where the Hamiltonian is zero, one can show that for virtually com- 
plete information destruction, Lindblad operators of dimension 2 log d are necessary 
and sufficient, where d is the size of the total Hilbert space (for simplicity we use a 
finite dimensional Hilbert space) We shall not reproduce the proof here, but it is an 
adaptation of the proof in fl3~5f . An example would be to have Lindblad operators 
which are projectors onto integrals of a complete set of local projectors, and their 
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conjugates e.g. projectors onto values of <fr n (x)7r m (x). If we demand that the map be 
only randomizing and not necessarily completely randomizing, then just over logd 
projectors are needed but they should be chosen at random. For Hamiltonians which 
are sufficiently strong or mixing, one can use less Lindblad operators for achieving 
information destruction. This is because in the interaction picture, the Hamiltonian 
evolution effectively acts a different Lindblad operator at each sufficiently large time 
interval. Understanding the time scale over which a particular Hamiltonian and set 
of Lindblad operators can achieve information destruction is a potentially interesting 
problem for future research. 

(5) Cosmological implications It is usually said that correlations in the 
cosmic microwave background (CMB) indicate that distant parts of the universe 
must have been in causal contact at some earlier epoch. Since the CMB is correlated 
over distances larger than the horizon, it is generally believed that the universe must 
have undergone some form of inflation, so that distant parts of the universe where in 
causal contact in the early universe. What has not been noted, is that correlations 
over distances which are space-like separated, do not necessarily imply acuasality. 
Indeed, as explained in Section [|, non-local signaling and non-local correlations are 
distinct, and the latter does not lead to a break down in causality. Correlations in the 
CMB need not indicate that a signal must have propagated between the correlated 
regions, or that the regions were ever in causal contact. One can have the creation 
of correlations in the CMB over space like distances in theories which do not allow 
signals to travel faster than light. Indeed, the theories we have examined generically 
create correlations over space-like separated regions, and we give a model theory in 
Appendix |B| which does exactly that. Even more intriguing, the correlations are 
created over all length scales equally. While this effect is usually small, it would 
be interesting to see whether it could be made to account for the observed CMB 
spectrum, in a model without inflation. 



11. Stochastic collapse models 

In this paper, our primary concern has been in decoherence models and information 
destruction. However, the exact same considerations also apply to spontaneous col- 



lapse models|36], [37], |38], |4(|. Such models are attempts to explain the collapse 
of the wavefunction, not epistemically, but as an actual dynamical process. It is not 
our intention here to advocate for such models. Rather, we merely wish to point out 
that some of the difficulties they suffer from are of a similar nature to the problems 
which plague decoherence models, and these problems may be addressed using the 
exact same methods outlined in the previous sections. 

In the classic GRW model 0, |37|, there is a probability per unit time that each 
particle is hit, meaning that it has a probability of being localized to a Gaussian 
wave packet at a particular point in space. The more particles at a particular point, 
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the more likely it is to be hit, thus macroscopic states, with a greater number of 
particles, are more likely to be localized. This leads to a possible explanation of why 
large objects behave classically (in the sense that they have a well defined position). 
One of the problems with the theory, is that it doesn't obey conservation rules. Other 
problems with the theory include a difficulty in constructing a field-theoretic version, 
since the non-conservation of energy leads to an instability of the vacuum. There 
has also yet to be a Lorentz invariant model. Finally, as noted briefly in Section [| 
the theory violates causality, on the scale of the size of the Gaussian wave-packet to 
which one collapses to. 

These problems can be solved using our relational approach. We will briefly 
discuss another difficulty - the so-called problem of tails. We shall not go into great 
depth here - we will simply present various collapse theories and then show how 
to modify them to address the problems of conservation law violation and lack of a 
field-theoretic model. Each collapse theory can be thought of as a different way to 
unravel the Lindblad equation. I.e. one has a theory which describes not just the 
density matrix, but each realization of the evolution in terms of a pure state which 
remains pure throughout the evolution. 



In the GRW [pq, 37] model, if a hit occurs on the i'th particle at point x, the 



wave function is multiplied by a Gaussian function 

G(q h x) = Kexp ("^fe - *) 2 ) (H-l) 

where a is some localization size, which can be taken a w 10 -15 cm, and qi is the 
position of the i'th particle. The probability Pi(x) that a particle is hit at point x 
is taken to be \{ipi(x)\ipi(x))\ 2 . The hitting occurs at randomly distributed times 
according to a Poisson distribution. 

We can conserve energy and momentum by not collapsing to location in terms of 
some absolute and unobservable external space-time, but rather in terms of a physical 
meaningful relational position. We also take the distribution in hitting times, to be 
given by a distribution in terms of some physical time r rather than an unobservable 
absolute time t. 

For example, consider the relational projector Qy as in Equation ( |3.1| ) acting on 
the i'th and j'th particle 

Qij = J dx\x(t)) ii {x(t)\®\x(t)} jj {x(t)\ (11.2) 

With probability \(ipij\Qij\ipi)\ 2 the i'th and j'th particle are hit, and well end up in 
the state given by Equation ( |11.2| ) appropriately normalized. 



To additionally conserve energy one can use the operator 

Qij = / dxdt\x(t)) u (x(t)\ ® |x(t))ii(s(t)| ® \T(t))(r(t)\ (11.3) 
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and the frequency of a particle being hit should be fixed in terms of the physical 
time r rather than the unobservable parameter t. 

The fact that the wave function has a probability of taking on values of x away 
from its peak is known as the problem of tails. I.e. the wave-function is hit by a 
sharply peaked Gaussian as in Equation flll.l ) and not a delta-function. The reason 



a Gaussian is chosen is that delta-function would lead to arbitrarily large deviations 
from energy and momentum conservation. This is particularly problematic in a field 
theory, where such an effect makes the vacuum unstable. 

Using relational observables may make the problem of tails less severe, since 
one is able to conserve momentum and energy even if the wavef unction is hit by a 
delta function. In such a case, one would still have large transfers of momentum 
from one particle to another, so one may still run into difficulties with experimental 
constraints, although often particles that are close together are held together by some 
potential. One can mitigate the energy and momentum transfer in a field theory, by 
using sufficiently gentle projection operators, as described in Section [|. The problem 
of tails is most sever in field theories, since collapse to a delta-function would make 
the vacuum unstable. As we saw in Section |] and [7], in the field theory, dynamical 
collapse models were stable in a field theory which has non-Gaussian collapses. 

If we still want to collapse the particles to Gaussians, then we could hit the state 
with 

G(q u qj) = Kexp{- — {q l - qj) 2 (11.4) 

which serves to localize particles, but around another particle which acts as a refer- 
ence frame, rather than around a point in absolute space. Such a theory is slightly 
non-local, with the degree of non-locality being on the scale of a. 

Let us now turn to another model, known as quantum state diffusion |39||. In 



such a theory, the evolution of the wavefunction is made to satisfy the Ito equation 

i#> = -iF^)dt+^x)( 2 ( M i) M i- M M-< M i) < Af i»iV'>^+x;(M i -(M j ))iv>^(t) 

3 j 

(11.5) 

where d£j(t) are independent complex and differential random variables representing 
a Wiener process. I.e. the mean of d£j(t), denoted by E(d£j(£)) satisfied E(d£j(t)) = 
while the mean of d£j(t)d£k{t) satisfies 'E[d£j(t)d£k{t)] = b~jkdt. This gives the 
evolution of a particular wavefunction, thus it is used for dynamical collapse models. 
When averaged over the random variable however, it simply gives the diagonal form 
of the Lindblad equation (|2.2|) . 

If we wish to have the Ito equation preserve conservation laws, we simply act 
as before, constructing relational operators Qi via Equation (|8.1|) and using them in 
place of the non-relational operators Mj. If we wish to conserve energy, we additional 
replace the parameter dt and t by a dynamical observable dr and r. For conservation 
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of energy and momentum, and Hermitian Q we would have 

|#) = -iH\^)dT-±J2@i ~ (Qj}) 2 \^}dr + Y,(Qj ~ (0;»|V>Ri(T) (11.6) 

o j 

Another dynamical collapse model is Pearle's continuous spontaneous localiza- 
tion model (CSL) which modifies the Schroedinger equation as 

^ME = -i#|V>,*>«-^ J dx[w(x,t)-2XA(x)] 2 }\^,t) (11.7) 

where A(^) could be any operator, but is usually taken to be the number of particles 
in a volume a m 10 _15 cm centered around x 

A{x) = J -^ Jl J dzN{z) ex P (-^^) (11.8) 

with N(z) the local number operator, w(x,t) a fluctuating field with probability 
density functional 

P T (w) = W (^ } T\^,T) W (11.9) 
and < t < T. The resulting wave function is then 

i r T 

|^,T) = Texp( / dxdt[w{x,t) - 2XA{x)] 2 ]) (11.10) 



4A 



T 



i.e. for a particular realization of the fluctuating field, it is hit by a random Gaussian 
centered around a particular value of A(x,t). 

An attempt at a relativistic model was made in 51]. As the author noted, it 
was unstable due to lack of energy conservation. We can make the theory stable 
and respect conservation laws if we replace in Equation Q11.7 ) the operator A(x,t) 
with our relational Q or Q, and t with the dynamical variable r to get a continuous 
localization model which respects conservation laws e.g. 

= -iHty, r) w - 7X>^) - 2A,Qi] V, r) (11.11) 

i 

Furthermore, since energy and momentum are conserved, one also may be free to 
have a sharper collapse 



d\ip,r) 



W^, t) w - ^9(w(t), Qi)W, r) (11.12) 

i 1 

where g leads to something other than a Gaussian (e.g. it could be the logarithm of 
a delta function). Furthermore, we saw in our discussion in Section |6|, that we can 
make the theory respect a form of Lorentz invariance. One approach is to make Aj a 
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dynamical variable which transforms as dt, the theory then respects minimal Lorentz 
invariance. 

Although the above model would solve the problem of tails, it is unclear whether 
it is within experimental constraints. Note that here, t) w could be a field wave- 
functional, thus one can have a field theoretic version of CSL. To minimize back 
reaction, we could proceed as in Section £| and take the local Lindblad operators 
Li(x) to be smeared local projectors e.g. 

Li(x) = J dlg t (l)Pi(x) (11.13) 

where the Pi(x) are projection operators onto values of some local field field or 
an operator acting on the local field (such as powers of the Hamiltonian) suitable 
smeared over a local test function. The gi(l) is a Gaussian or other distribution 
peaked around some value of the field. One then group-averages this operator to 
construct the Qi used in any of the collapse models (while some use only one Q, 
there is no obstacle for using many). 

We have not said much about what observable the Lj(x) should represent. Indeed 
there is not a good reason to a priori prefer any particular one. Traditionally they 
have been taken to be the positions of particles, but they could also be energy density, 
or some other such observable. For example something along the lines of integrals 
of the local number operator or energy density, and powers of it, as discussed in 
Appendix are tempting choices since the strength of decoherence in the Lindblad 
equation is proportional to the number of particles, and thus produces the desired 



effect that macroscopic objects are more likely to be effected (c.f. f43l). 

Finally, we note that although here we have discussed the set of relational ob- 
servables in the context of spontaneous collapse models, the discussion applies to 
other interpretations of quantum mechanics. For example, in Bohmian theory, the 
particle trajectories do not conserve energy or momentum, and similar considerations 
might be applied there in order to find the physically relevant trajectories which do 
conserve energy or momentum. Even in a many-worlds interpretation of quantum 
mechanics, one might consider our relational observables as describing the relevant 
observables in a particular branch that is observed. 



12. Conclusion 

We have seen that it is possible to construct a fundamental theory which allows for 
information destruction, while still being causal and respecting conservation laws. 
Whether such a theory is consistent with experiment, or is in fact realized is an open 
question. However, the mere possibility of such a theory is enough to justify future 
study, especially in light of the motivation coming from black hole evaporation. This 
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leaves many open questions, some of which we shall turn to after briefly summarizing 
our results. 

We followed a natural generalization of unitary theories to non-unitary theories, 
through the Lindblad equation of 



-£ = -i[H,p]-^^2 7ij (M]M i p + pM}M i -2M i pM}) . (12.1) 



dt 2 



This is the most general evolution one can have which is a completely positive map 
(doesn't change the meaning of local density matrices), and is Markovian. However, 
in order to ensure that conservation laws were respected, we replaced the Lindblad 
operators Mj with the relational operators 

Qt = J dgU(g)M^ 2 ' 3 - ) U^g) (12.2) 

where we average over some unitary representation of a group so that our evolution 
conserves the generator of that group. This allows us to respect the traditional con- 
servation laws. We also introduced extra fields 1,2,3... which allows our operator 
Qi to be nontrivial, as it can describe relations between different fields. We ar- 



gued in Section |8J] that such operators allow for the destruction of the information 
concerning the sorts of observables we make in the lab. 

Energy conservation needed to be treated slightly differently. The first step is the 
same as above, and involved group averaging, using the Hamiltonian as a generator 
of time translations. The resulting Lindblad operators conserved energy but were no 
longer local. In order to make them local, we imposed time-translation invariance 
on the state, 

d -l = (12.3) 
dt y 1 

which yielded an equation akin to the Wheeler deWitt equation. This equation was 

local and conserved energy, and was non-trivial in the sense that one could describe 

evolution in terms of correlations between observables and a clock. 

In Section |5|, we showed that this evolution respected causality. This was done by 

taking our relational Lindblad operators to be integrals of local operators, and taken 

them to be either Hermitian, or to come in a pair with their Hermitian conjugate. We 

also distinguished another effect - the creation or destruction of spatially separated 

correlations. This an effect was not seen in the non-relativistic limit, but was seen 

for relativistic fields. Fortunately, it does not result in a break-down of causality. 

Certainly this evolution of correlations is non-intuitive, but it cannot be used to 

signal superluminally. We discuss the effect of correlations in a model given in 

Appendix [FJ. 

In Section |9|, we also took on the question of to what extent any symmetry 
should imply a conservation law. We derived a generalization of Noether's theorem 
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and the continuity equation for Markovian non-unitary evolution. This placed a 
restriction on the Lindblad operators - they either had to be Hermitian (in which 
case, the symmetry implied the conservation law), or they had to be creation or an- 
nihilation operators. Combined with the locality constraints, this appears to restrict 
the violation of conservation laws to be mild. It remains an open question whether 
further restrictions exist which would require exact conservation from the symmetry, 
or whether such a restriction ought to be imposed to keep the theory finite. 

In Section || we discussed the fact that our theories were invariant under the 
orthochronous Lorentz group, but proved that it was impossible to get Lindblad 
evolution where the dissipater term transforms like a four-vector. We also show in 
Appendix y, that the theory could be made time- symmetric. 

Finally, in Section [II], we discussed application of the relational approach taken 
here, to stochastic collapse models believed by some to describe the emergence of 
classicality from quantum mechanics. While in Section [10] we discussed toy models 
for coupling the theory to gravity. 

We have left open many theoretical questions, and we take the opportunity in 
this conclusion to summarize some of them. From an experimental point of view, 
one might test some of the models, by looking for momentum exchange between 
pairs of particles. Such an effect may lead to spontaneous ionization of atoms, as an 



electron and nucleus are given a kick in the opposite direction [ 44 1. The existence of 



such an effect, and it's magnitude, depends on which model is chosen. In some mod- 
els, the effect can be made very small making it unlikely that current experimental 
constraints could rule out the theories considered here, although they may restrict 
the parameters in some models. Likewise, long range creation of correlations, or 
abnormally fast decay of correlations is a signature of all the models, but we suspect 
its effect would be very difficult to measure. 

From a theoretical point of view, some of the more pressing open questions con- 
cern the richness of local relational observables, as noted in Section |8.1] . To what 
extent can such observables describe the world around us, or destroy the required 
information in a black hole. To explore this, it may be interesting to look at numer- 
ical simulations of some of the models discussed in Sections ^O] and Appendix [A], 
particularly those which result in effective potentials which depend on distance. 

The constraints due to causality discussed in Section [| and Appendix [D] were not 
fully explored, and more work is needed in determining which Lindblad operators 
are consistent with causality, and lead to renormalisable theories. Even constructing 
additional models to the one we described would be useful. 

In terms of coupling the theory to gravity, we only mentioned possible toy models, 
and it would be interesting to pursue this further. Especially with respect to theories 
which don't conserve the energy-momentum tensor, since the difficulty in coupling 
such theories to gravity might preclude such theories. With regard to gravity, we 
noted that correlations due to fundamental decoherence could be imprinted in the 
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cosmic microwave background, and it would be interesting to explore this further, 
possibly as an alternative to inflation. 

With regard to the generalization of Noether's theorem proven in Section ^ it 
would be interesting to see whether it could be applied to any condensed matter 
systems. Nothing about it depends on whether the non-unitarity is fundamental, we 
just require that some symmetry be respected, perhaps only approximately. 
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A. Models 

In the body of this paper we have tried to keep our considerations as general as 
possible, using particular models only for illustrative purpose. In this appendix we 
will briefly sketch a few particular models, which we have either already used, or 
which we have found useful. We will not explore them in any detail, but rather hope 
that they might prove useful for future research. 

We will first introduce Lindblad operators based on a local number operator. The 
motivation for this is that it is a field theoretic model, whose non-relativistic limit 
gives the quantum mechanical model of Equations ( |3.1[ - |3.2|) . We wish to construct 



a field theoretic model which has as its limit, the quantum mechanics model of 
Equations fl3TT| - |3.2|) . To this end, consider a scalar field and define 



vb(x) = <t>{x) + in(x) (A.l) 

In the non-relativistic limit ij}'(x) acts as a localized creation operator i.e. it creates 
a single particle which is localized around the point x 

|l.>=^(x)|0> (A.2) 

We wish to construct a field theoretic model which has as its limit, the quantum 
mechanics model of Equations ( |3.1| - |3.2| ). To this end, consider a massive scalar field 



and define 

vj(x) = <f>(x) + iir(x) (A. 3) 
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The operator ip<(x) acting on the vacuum creates a single particle state 

\l x )=if)\x)\0) (A.4) 

which is localized around the point x. 

For modes with k <C m we get if>(x) oc J dke lkx a k and hence for non-relativistic 
modes if> and ^ acts as a non-relativisitic annihilation operators. 

One can verify that \l x ) transforms correctly under translations i.e. 

Ui\l x )(l x \U} = \l x+l )(l x+ i\ (A.5) 

where C/j = exp(iZP) and P is the field momentum operator. 
Furthermore, one also has 

(1*11*+,) = f dk(2 + u k + —)e ikl -> (A.6) 

27T ,/ L) k 

for I »> 1/m. In this case there is a region for which u k + — rs constant and the 
integral tends to zero. Hence the state indeed describes a single particle state 
that is localized around x over a region of size 1/m. 

We can then define a number operator as N(x) = ip j; (x)if>(x) with if) defined as 
in Eq. (|A.3|) . Now imagine two fields, so that we have Nx(x) and N 2 (x). Then take 

Q = j ' dxN^x) ® N 2 (x) (A.7) 

as the Lindblad operator we are interested in. In the single particle non-relativistic 
limit, i.e. the probability of getting N±(x) > 1 or N 2 (x) > 1 is very small, this 
Lindblad operator clearly reduces to that of Equation (|3.1| ) 

Q = J dx\x){x\ ® \x){x\ 

We thus see that we recover the quantum mechanical model in the non-relativistic 
limit. 

If we wish to have a field act as a reference frame, then another method is to 
use coherent states. In this limit, we can have a fully field theoretic model which 
can reproduce the models of Equation ( |3.3| ) used to show that one could decohere a 
field at a point and still conserve momentum. The coherent states behave in some 
sense like a classical system, and thus behave like localized particles. Recall that for 
a harmonic oscillator we can construct a coherent state 

\g) :=e 9a+9 * at |0)^o (A.8) 
from \0)ho t ne ground state of the Harmonic oscillator such that a\g) = g\g), (g\aJ = 

g*{g\- 



49 



Since each field mode is a harmonic oscillator we can construct something like 
\^coh) = ®k\gk(%o)) for the field 

\^ coh ) = e /*(fl(*)«k+fl*(*)«t)|o) (A.9) 
If we compute the expectation value of the field at a point it will be 

/rlk 
-—(g(k)e^ + g*(k)e~^) (A.10) 

so now cifc is an eigenstate and the expectation value can be chosen by taking suitable 
g(k,x ) to be concentrated around x , so that {</>(x)) = g(x — x ). 
Suppose we have two fields, and choose 

Q = J dxM*)Mx) (A.11) 

although any local operator L(x) can be used in place of 2 {%)■ Let us look at the 
dissipater term 

V{p) = -i\[QAQ,p\ (A-12) 

if we assume that p ~ p\ ® p 2 and take e.g. p 2 = \^ C oh) ($coh\ than under the 
assumption that the evolution still keeps p 2 close to the chosen coherent state the 
dissipater term is roughly give by 

V(p) « -7^ J dxdx'g{x - x )g(x' - x o )[0i(x)0i(x'), p]] (A.13) 

which is exponentially small if x and x' are not close to Xq. 

If it is possible to take the limit that g(x — Xq) are delta functions we get that 
effectively the dissipater term is 

^(p)^-^7[0 2 M,[0 2 (x o ),p]] (A.14) 

a decoherence term at a point, and yet we did not violate momentum conservation. 

Turning back to models which use number operators such as those of Equation 
QA.7 ), we might generalize this to consider higher and higher powers of N(x), i.e. 



Qmn = J dxN™{x) ® N^{x) (A. 15) 

in order to get at the full richness of the theory. 

In such a case, one can even consider only one field, since such a theory can be 
made self-referential i.e. if we have two particles and only one species, then the set 
of Lindblad operators 

Q m = I dxN m {x) (A.16) 
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will distinguish between the cases where the particles are coincident or not. In 
the limit where we only have two localized particles, one only needs Q2, since non- 
coincidence gives Q2 = 2, while coincidence gives Q2 = 4. In effect, one has a know 
background field O which acts as a reference frame, and one treats the remainder, 
5(fi = — 4> as a perturbation which is being measured or decohered. Thus to first 
order 

Q 2 « N a + 2 / dxN {x)5N{x) (A. 17) 



which, when N is fixed, acts very similar to the two field model described by Equa- 
tion ( |OD 

It is also very tempting to consider models which have Lindblad operators as 
integrals of powers of the Hamiltonian 

Q m = J H m (x)dxdt (A.18) 

However, it remains to show that such a theory is renormalisable. Nonetheless, 
such models may be particularly interesting for stochastic collapse theories, since the 
greater the energy density, the greater the probability of collapse. 

To avoid back-reaction, and keep the theory finite, we might want to smear over 
field observables. This can be done both over a small region of space, and also over 
values of the field. We discussed this at the end of Section (|, and so shall not repeat 
that discussion here. 



In Section |8J] we discussed why one might want to consider models which have 



tiny violations of non-locality such as those given by Equation 

Qiap = J dxA a (x)B{3(x + I) 

we might then integrate over I to create a smeared Lindblad operator over some small 
range of I. 

An analogous construction can be considered as follows. Let E x = \l x )(l x \ be 
the projection-like operator on a localized single particle defined above, and E x be 
the same type of operator for a second field. Now define the Lindblad operator 

Q = J dxdx'E x E x ,h(x-x) (A.19) 

where h(x — x') is a "hat" shaped function, h — 1 for \x — x'\ < A. 

Now we have two scales 1/m the "size" of the localization of E x , and A the size 
of the smearing. An interesting case is A ^ 1/m, where we smear over scales larger 
then the Compton wave length of the two fields (i.e. we require A >> max(l/m, lm)) 

Under this assumption Q 2 w Q. I.e. in the limit it becomes a projector. 

Q 2 = J dxdx'dydy'E x E y E x/ E y ,h(x-x')h(y-y') (A.20) 
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but since most of the integration is over regimes for which x — x' ^> 1/m and 
t/-!/'>l/mwe can replace 



E x E y « 5{x-y)E x 



(A.21) 



E x iE y i « <J(x' - y')^ 



(A.22) 



so we get 



Q 2 



dxdyE x E x ih 2 (x — x') 



(A.23) 



but since h — 1 we obtained Q 2 = Q (equality of course only in the limit of mA — > oo, 
but we can perhaps expect, that Q 2 tends to Q up to a correction that dies of like 
exp — (mA)) 

The key point is that due to the smearing function the field which acts as a 
reference frame is "too big" to resolve close to the Compton length scale. Once we 
do that, we retrieve that Q is a good projector even in the fully relativistic case. 
This smeared projector commutes with Ptotai since the function h(x — y) depends on 
relative coordinates. 

On a completely different note, recall in Section |8.1| that there are models in 
which one can add potentials to discriminate between particles, or localized fields at 
different distances. I.e. the Hamiltonian terms might become very important when 
one has an effective potentials V(x — y) which depend on the distances between 
objects and couples to other degrees of freedom (such as local spins). This can serve 
to distinguish distances between particles. Since analyzing such models probably 
requires numerical calculations, we will not treat such models in more detail here. 

Finally, another set of possibilities, is to take operators which are coherently 
group-averaged as an example, first recall that our coincidence Lindblad operator of 
Equation ( |3.1| ) was 



which is an incoherent mixture of coincidences, representing the fact that one lacks 
knowledge of the absolute reference frame x. However, one could consider the func- 
tion generated by coherent group-averaging 



Just as before, this commutes with P (it is an eigenstate of P, rather than a mixture 
of eigenstates. The projectors Qk onto \ipk) could then be used in the Lindblad 
equation. These operators are rather trivial, but this appears to be a function of the 
simplicity of the case being considered. 



Q 



dx\x)(x\ <S> \x){x 




(A.24) 
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For the field, one can consider something similar, by looking at operators build 
out of a coherently group-averaged operator 

M t = J dgdg'U{g)M l U\g') (A.25) 

This operator is invariant under incoherent group averaging, and will thus conserve 
the required quantity. However, it is not always local. 



B. Correlations in a particular model 

Here we examine a particular model and show that the creation or destruction of 
correlations occurs with very small probability on typical states. I.e. it takes place 
on a vanishingly small part of the full state space. For simplicity, we will consider a 
quantum mechanical version of the field theory operator of Equation (P-l| ), however, 
our considerations will apply to a field theory. To this end, let us take Pi(x) to be 
a set of smeared projectors (one might be inclined to imagine that they are akin to 
projections onto the values of some field observable <p{x)) in a discretization of <p(x)). 
Let us take as an example, a model in which we can write 

Qi= [ dzPfc) (B.l) 



with Pi(z) some projector or smeared projector. Next, consider only two points in 
space, x and y, and define 



>'!) 



P(x) ® I(y) + I(x) <g> Pi(y) . (B.2) 



Then dropping the Hamiltonian term for simplicity we get that evolution of an 
operator A(x) is local, and of Lindblad form. I.e. 

j t A(x) = -^[Piix), [P{x), A{x)}\ (B.3) 

as we expect from Section |. However, the correlation A(x) <S> B{y) of two distant 
observables behave as 



j t [A{x)®B{y)\ = -±Y,li[QiAQiM{x)®B 

i 

= -\j2j dzdJ'y i [P i (z),P i (z'), [A(x) ® B 



2 

i 

-A{x) ® B(y) + A{x) ® -B(y) + V(A(x) ® B{y)) (B.5) 
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where V is the term which violates the product evolution form of Equation ( 5.21 
and is equal to 

V=- J2l p i(xl A{x)\ ® [ p i(vl B (V)} (B.6) 

i 

this violating term looks a bit like a Lindblad form, as it can be written as 

v = - J2l p i(xl My), M*) ® B (y))} (B-7) 

although the matrix of projectors is not positive. 

This terms is zero on the subspace where (p(x) ^ f(y)- he. consider the projector 
Pdis onto this subspace 

p*= E p i 1] ( x ) p j 2) ^) p i 1] (y) p i 2 \y) (b.8) 

then 

l/(P dis A(or) ®B(y)P ( «.) = (B.9) 

For a finite dimensional Hilbert space of dimension n at each point in space, the 
subspace in which V is non-zero grows as n while the total space grows as n 2 . Thus 
as n goes to infinity, the fraction of the Hilbert space in which one finds creation 
or destruction of correlations goes to zero. Thus for typical states, destruction or 
creation of correlations will be a very small effect. If we increase the number of fields 
and how fine grained the projectors are, this term becomes rarer still. 

It is worthwhile to see the effect of the term V(pab) on the evolution of a state 



and how it violates the product evolution form of Equation ( 5.21|) . Let us consider a 
finite dimensional system and density matrix at points x and y decomposed in terms 
of the basis \i x )(i x \ = P%(x), 

P = <7ij,ki\ix){jx\ ® \ky){ly\ (B.10) 

The solution of Equation ( |B.4|) with and without the product violating term V is 
given in Table |B| for 7$ = 7. Note that the evolution is only non-product for terms 
\ix)(jx\ <S> \iy){j y \- These off-diagonal elements decay faster than they would under 
purely product decoherence which obeys Equation ( |5.21| ). 

This shows how correlations can decay faster than one might otherwise expect. 
One also finds creation of correlations. As an example, take as an initial state that of 
Equation ( |B.10|) and with all cr^^i equal. This is an initial uncorrelated pure state, 



with each party's state in a superposition over all basis states, i.e. we will take the 
state at x to be = \i) / y/n and similarly for the state \iPb) at y. Now, we 
know from the local evolution law of Equation ( |B.3|) that all off-diagonal terms of the 
local states will decay, and the local states will be maximally mixed. I.e. the state 
at x will evolve to pa = Yl K)(^l/ n > an d similarly for p# at y. The local entropy of 
each state is thus logn. However, from Equation ( B.2 ), we see that there are only n 
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condition on i,j,j, I 


with V(a) 


without V(a) 




&ii,kk{t) = Oii,kk{fy 


&ii,kk\t) — cr ii,fefc(0) 


k I 


cru,ki(t) = e-^a^kiiO) 


cru,ki{t) = e-^a^kiiO) 




crij,kk(t) = e~^a ijt kk(0) 


(?ij,kk{t) = e" 7 V iiifcfc (0) 


k ^ I, j,i ^ k, j ^ I 


VijMit) = z~ 2lt Vij,ki(®) 


crij,ki{t) = e'^VijM®) 









Table 1: Evolution of correlations of the density matrix p = (Tij t ki\ix){jx\ ® \k y ){l y \ with 
and without the product violating term V(a) 



projectors P® y . Thus the total state will only be decohered into n states, rather than 
n 2 possible states. In particular, although each local density matrix is maximally 
mixed, the total density matrix is decohered into the states 



\il> i ) = y %2\ij + ki)/Vtoi (B.ll) 

jk 

which can be seen from the fact that these are the states which survive when we 
apply the Lindblad operator to both sides of the initial state 

^MbKMbIpt 1 = mm (B.12) 

To put it another way, each Lindblad operator P® y is degenerate, selecting states 
where one of the two sites is in state but the other site can be in any state. Thus 
terms which are superpositions over these degenerate states will survive (the states 
of Equation ( [B.ll )). Thus although the local states will look maximally mixed, 



the total state will contain coherences. As a result, we have that the entropy of 
each system is n while the entropy of the total system is also n. Since the mutual 
information I(A : B) between two systems A and B is 

I(A : B) = S(A) + S(B) - S(AB) (B.13) 

where S(A) and S(B) is the Von- Neumann entropy S = —trplogp of the state at 
x and y respectively and S(AB) is the entropy of the joint system. In the present 
case, the entropy of the local states is log n, while the entropy of the total state is no 
larger than logn (the number of states the system is decohered into), and thus the 
mutual information is at least log n and approaches it for large n since then is when 
the states of Equation (|B.11| ) become orthogonal. 



C. A time-symmetric formulation 

A Lindblad equation has an implicit boundary condition at t — 0, because if you 
evolve the equation backwards in time past t = the equations no longer correspond 
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to a CPT map i.e. the density matrix will evolve into something which is not a 
density matrix. This is reasonable, since at t = we assume some initial condition 
(such as the state being pure), and then evolve forward in time. 

As a result, it is often said that the Lindblad equation gives an arrow of time, 
since as we evolve it, a state become more mixed. I.e. it is said that since the evo- 
lution gives an entropy increase, it must necessarily give an arrow of time. However, 
this simply reflects a choice of time-asymmetric boundary conditions - i.e. we start 
with a pure state at t Q and find that we don't have complete knowledge of its state 
at future times. Yet, given a pure state at t we could just as well have tried to 
retrodict what the state would be at earlier times. In this case, if the fundamental 
evolution is not unitary, then we will also not have complete knowledge of the state 
in the past. Such a situation has been contemplated for retrodicting the results of 
past measurements from future measurements, both in standard quantum theory [[45 



and in open quantum systems p6 



Here, we are interested in constructing a full retrodicting Lindblad equation for 
the density matrix itself. This has an entropy increase backwards in time, which is 
what one expects since for an initially pure state, you can't predict with certainty 
which state led to the current pure state. Thus we see that entropy not only increases 
into the future, it also increases into the past - therefore, it does not give an arrow 
of time. 

The retrodiction equation can be constructed from the forward evolving equation 
by changing the sign of the dissipater T>(p) so that if we were to take dt — > —dt, we 
would get Lindblad evolution backwards in time which is valid for earlier times. 

| = -^p]+^£ Hj(L}L iP + P L)U - 2L lP L)) (C.l) 

ij 

We could also derive the above equation by considering a microscopic derivation, 
where we have a unitary interaction and trace out an environment. I.e. starting 
from the usual forward evolving equation 

= -i[H E + H S + H ES , P se] (C.2) 

which leads to the usual forward evolving Lindblad equation, we take the backward 
evolving unitary evolution 

^ = i[H B + H s + H ES , Pse] (C.3) 
to get the backward evolving version of Equation ( |C.1| ) 

^ = i[H, p]-\Y, iniA^P + P L ) Li ~ 2L ^ ( c - 4 ) 

ij 
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where we have removed the labels S and E for system and environment. In these last 
two equations, a positive dt represents a step back in time, so we then take dt — > —dt 
to ensure that time is defined in the standard way, thus arriving at Equation ( |C.1| ). 
We can combine the integral version of both the forward and backward evolving 
Lindblad equation to define the evolution over all times t given some state at t = 



Jo \+S a dt\Y.ML\LiP + pL\L i 



D. No instantaneous signalling and causality 

In Section [|, we saw that as long as we chose the Lindblad operators to be integrals of 
local operators and either Hermitian, or coming in pairs with its Hermitian conjugate, 



then the evolution of of Equation (4.2) will be causal. In the present section, we note 
a more general locality condition for the Lindblad operators which is exact, but 
under the weaker requirement that the equations of motions of local observables are 
local at an instant, i.e. their evolution only depends on the value of local fields. 
This is equivalent to saying that there is no instantaneous signalling. I.e. for a local 
Hermitian observables A(x) 

= f{x) (D.l) 

with f(x) some local operator. This need not guarantee that signals cannot propagate 
at speeds faster than light, but it is a necessary condition. 

We can also write this condition as a vanishing of the equal time commutator at 
distant locations: 

[ ? t , B(x', t)} = C(x, t)5{x - x') (D.2) 

For dA(x,t)/dt we can substitute the right hand side of Equation ( |4.2| ), and expand 
the Lindblad operators as 

Qi= [ dxLi(x) (D.3) 



where we make no assumption about Li(x). This gives the general condition that 

for X ytz x' 

j dzdy[L\(y)[U(z),A(x)} + [A(x), L\(z)]L t (y), B(x')] = (D.4) 

It comes as no surprise, given the results of the previous Section [5], that one way to 
satisfy this condition is to require that for any local Hermitian operators A(x) and 
B(x): 
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1. (a) [Li(z),A{x)] = Di, A {x)6(x-z) 
(b) [D i>A {x),B{z)](x5{x-z) 

and that either 

2. Li (x) = L\{x) or 

3. for any ^f~flL ifl (x) = iff^Ll^x), 

where we have used a double index for the last condition in order to ensure that 
the Lindblad operators come in pairs with their Hermitian conjugate as in Equation 
( |5.18|) . Imposing condition 1 (a) and (b) (i.e. that Lj(x) is a local operator) on 
Equation (|D.4|) gives 

£ 7 < [D[ B {x')D hA {x) - D hB {x')D[ A {x)) (D.5) 

I.e. that li^l B {x')D itA (x) is Hermitian for all A(x), B(x). This can be satis- 
fied by taking the D^ A (x) anti-Hermitian or Hermitian, which is equivalent to taking 
the Li(x) Hermitian (Condition 2). It can also be satisfied by instead taking Condi- 
tion 3. 

For the more general condition of Equation (p.4f) , it is not clear what other 
ways exist where one could satisfy it nor whether other conditions would violate full 
causality. We conjecture that Conditions 1 and either 2 or 3 are required. 
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